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EDITORIAL 


À necessidade e o desejo de tornar acessível, aos profissionais e estudantes de enge- 
nharia, o conhecimento de novas técnicas e contribuições teóricas importantes resultantes 
da investigação realizada tanto a nível nacional como internacional manifestou-se há vários 
decénios nas escolas de engenharia, 


Assim, em Outubro de 1915 aparece em Lisboa uma revista designada por TÉCNICA 
INDUSTRIAL. Face à instabilidade política então reinante esta revista teve uma vida bas- 
tante efémera traduzindo-se apenas na publicação de pouco mais de uma dezena de nú- 
meros. 


É só em 1925 que se consegue a publicação de uma nova revista que, em princípio, 
deveria vir dar continuidade à ideia presente à fundação da TÉCNICA INDUSTRIAL. 


Dessa nova publicação foi posto à luz do dia o primeiro número em Dezembro de 
1925 ostentando uma nova designação — TÉCNICA—, 


Desde então ela vem sendo publicada regularmente tendo apoiado e divulgado muitos 
e importantes trabalhos especialmente no campo da engenharia. No entanto esta regulari- 
dada nem sempre se pôde considerar óptima o que aconteceu por motivos vários de todos 
já conhecidos e várias vezes explicitados. 


Merecem, apesar de tudo, referência os encerramentos das instalações a que as auto- 
ridades por diversas vezes nos impuzeram. Isto originou forçosamente uma grande descoor- 
denação quer interna quer a nível de contactos externos o que só foi possível ultrapassar 
com grandes esforços não só dos colaboradores mas também dos autores. 


Para avaliar da gravidade destes acontecimentos bastará notar que em cinquenta 
anos de actividade, que deveriam corresponder à publicação de quinhentos números, só se 
nos tornou possível a publicação de quatrocentos e trinta e um como se pode observar 
pelo número da presente revista, 


A ambos queremos inequivocamente expressar o nosso reconhecimento na certeza de 
que os seus esforços não foram em vão, 


No momento em que comemoramos o nosso primeiro CINQUENTENÁRIO não nos 
poderiamos furtar a esta simples homenagem fazendo votos que doravante as nossas 
actividades possam prosseguir sem graves interferências permitindo o desenvolvimento 
normal da tarefa em que nos empenhamos. 


A Técnica 
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O artigo DO CONHECIMENTO E DA ACÇÃO de ANTÔNIO GOUVÊA PORTELA foi 
publicado inadvertidamente com alguns erros tipográficos, 


Ao autor e leitores apresentamos as nossas desculpas. 
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Oscilações auto-excitadas num circuito 


condensador - termistor (') 


RESUMO 


Faz-se a análise do comportamento dinâmico da 
associação em paralelo dum condensador com um ter- 
mistor alimentada por um gerador ideal de corrente, e 
em particular da possibilidade de ocorrência de oscila- 
ções auto-excitadas. Comparam-se resultados experi- 
mentais com os resultados da simulação numérica do 
modelo teórico proposto. São descritas as experiências 
que permitem determinar os parâmetros do termistor 
8 uma montagem destinada a simular os elevados valores 
de capacidade que foram necessários nas experiências. 
Utilizando a simulação numérica, faz-se o estudo e in- 
terpretação fisica do comportamento do sistema para 
diferentes valores da capacidade do condensador. Apre- 
senta-se, para caracterização do termistor no dominio 
incremental, um esquema eléctrico equivalente deduzido 
das equações diferenciais não-lineares que regem o 
seu comportamento. 


1 — INTRODUÇÃO 


O termistor é um elemento óhmico com um par de 
terminais acessíveis, cuja resistência é fortemente de- 
pendento da temperatura, descrescendo com o aumento 
desta. 

Em meio ambiente a temperatura constante, o auto- 
-aquecimento produzido pela passagem da corrente con- 


J. F, Borges da Silva 
Doutor-Engenheiro |. S.T. 

Prof. Extraordinário do |. 5.T. 
J. M. Viegas Gonçalves 

Eng. Electrotécnico 1. S.7T. 
Assistente do |.5.T. 


ABSTRACT 


Tha dynamic behaviour of a condenser and a ther- 
mistor connected in parallel with an ideal current source 
is analysed, particular attention being given to the possi- 
bility of self-excited oscillations. Results from experi- 
ments are compared with results obtained by numerical 
simulation of the theoretical model proposed. Experi- 
ments required for the determination of thermistor para- 
meters and a method for the simulation of large capaci- 
tance values are described. System behaviour was further 
investigated for various capacitance values using numeri- 
cal simulation, the results being interpreted physically. 
A thermistor small signal equivalent circuit, derived 
from the non-linear differential equations describing its 
behaviour is presented, 


duz a uma característica estática corrente-tensão com 
um troço decrescente, Fig. 1. Dependendo a tensão U 
dos valores da corrente | duma forma univoca, um ponto 
de funcionamento na região descendente da caracteris- 
tica será estável, se a resistência do circuito exterior for 
superior à resistência incremental do termistor (em mó- 
dulo), no ponto de funcionamento, isto é, se o regime 
for do tipo «corrente imposta». No troço descendente 


(1) Publicação do Laboratório de Medidas Eléctricas do 1.5. T. integrada no Projecto de Investigação TLE-5 subsidiado pelo Ins- 


tituto de Alta Cultura. 
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da característica o termistor é, todavia, potencialmente 
instável com propriedades, sob este ponto de vista, aná- 
logas às do arco eléctrico. 


Em regime dinâmico o tempo intervém necessaria- 
mento na relação corrente-tensão do termistor, pois que 
em virtuda da inércia térmica, a temperatura que con- 
diciona o valor da resistência, não segue instantanea- 
menta as variações da potência de auto-aquecimento 
P = UI. Verifica-se que este efeito dinâmico se pode 
traduzir, em regime incremental e do ponto de vista 
eléctrico, por uma indutância equivalente, 


Os factos anteriormente apontados sugerem a pos- 
sibilidade de se obterem oscilações auto-excitadas, asso- 
ciando um condensador em paralelo com o termistor, 
Fig. 2, e alimentando o conjunto por meio dum gerador 
ideal de corrente contínua, que imponha um ponto de 
funcionamento em repouso no troço decrescente da ca- 
racterística estática. 


Fig. 1: (a) característica estática dum termistor. 
(b) Característica externa do gerador (caso estável). 


Verifica-se efectivamente que, nas condições enun- 
ciadas, o ponto de funcionamento em repouso se torna 
instável para valores da capacidade superiores a um 
certo valor crítico, dando origem espontaneamente a 
oscilações de amplitude crescente, que vêm a estabilizar 
num regime periódico de amplitude constante. 


Em regimes próximos do crítico as oscilações apre- 
sentam carácter quase sinusoidal, degenerando em os- 
cilações do tipo «relaxação» para regimes fortemente 
sobre-críticos. Estes últimos são, quanto às formas de 
curva, semelhantes aos que se produziriam substituindo 
no esquema da Fig. 2 o termistor por um tubo de gás. 

A ordem de grandeza das constantes de tempo 


inerentes ao processo de variação da temperatura dum 
termistor, mesmo de pequenas dimensões (100 mm), 
conduzem a oscilações com periodos na casa da dezena 
da minutos, 

Ocorrem na literatura referências a algumas das 
matérias abordadas no presente trabalho, de que se 
citam [1]. [2]. 


Fig. 2: Circuito condensador-termistor Tr- Termistor. 


2 — OBJECTIVOS DO PRESENTE TRABALHO 


1) Estabelecer por análises do esquema da Fig. 2 
o sistema de equações diferenciais descritivo do seu 
comportamento dinâmico e investigar as condições de 
instabilidade, 


2) Obter por via experimental os parâmetros rele- 
vantes para a descrição do termistor. 


3) Simular numericamente o comportamento do 
sistema da Fig. 2 com base nas equações de 1) e va- 
lores obtidos em 2), pelo emprego dum algoritmo de 
integração de equações diferenciais não lineares. 


4) Confirmar a validade da simulação, comparando 
os seus resultados com medidas efectuadas sobre um 
sistema físico com a configuração da Fig. 2. 

5) Utilizar a simulação numérica para investigar 
outras situações, para além daquelas que foi possível 
observar directamenta nas experiências. 


3 — ESTABELECIMENTO DAS EQUAÇÕES DINÂMICAS 
3,1 Descrição do termistor 


A substância dum termistor pode ser concebida 
como um semicondutor com impurezas." A variação da 
condutividade com a temperatura traduz essencialmente 
a variação da densidade de portadores de carga liber- 
tados por ionização dos centros de impureza, À lei que 
defino a dependência da condutância dum termistor re- 
lativamente à temperatura apresenta assim a forma expo- 
nencial típica dum fenómeno de ionização, [3]. 


B 
G(T) = G, exp = 
dl ci E (3.1) 


(1) Exemplo típico são os termistores fabricados à base de óxidos dos metais do grupo do ferro, como (0--), ([Fe+++ |). 


consistindo as impurezas na substituição de alguns átomos de Fe por Ti, formando centros (07- dh, Fa +++ Ti+t+tT; om que o 
electrão transferido de Ti para Fe está fracamente ligado. A temperatura ambienta a libertação destes electrões origina uma condução 


do tipo n, 13). 
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onde: 


T — Temperatura absoluta em “K. 

B — Constante em “EK. 

G -— Condutância do termistor em 5. 
Go 


— Constante em 5. 


Pela especificação das duas constantes B e G, 
fica completamente definido este aspecto das proprieda- 
des do termistor. 

Do ponto de vista térmico, hã a considerar a quan- 
tidade de energia libertada por efeito de Joule na uni- 
dade de tempo, 


Ul = U"G(T) (3.2) 


o a função que traduz a transmissão de calor ao meio 
ambiente, (ar tranquilo), que se admite a temperatura 
T, constante. 

Considera-se a transmissão D como uma função da 
diferença de temperatura entre O termistor e o meio am- 
biente com a forma, 


DATT, )=S. (TT, )+HO (TT, )º (3.3) 


onde: 

D — Ritmo da transmissão de calor em W. 

S — Constante em WKT| 

O — Constante em WºK”* 

B transmissão de calor por condução pura, daria 
origem apenas ao termo linear de (3.3). Verificou-se 
experimentalmente a necessidade de juntar um termo 
quadrático de correcção, o que evidencia a influência de 
outros mecanismos de transmissão de calor, por exemplo 
a radiação, com importância para as diferenças de tem- 
peratura mais elevadas. 

Em face das relações (3.2) e (3.3), o balanço tér- 


mico instantâneo do termistor pode ser expresso pela 
equação, 


dT 
k—> =—D (TT, ) + U'G(T) 


dt (3.4) 


onde: 


t — Variável tempo em s, 


k — Capacidade térmica do termistor em JK”! 


3.2 Descrição do circuito eléctrico 
Reportando-nos ao esquema da Fig. 2, obtém-se 


imediatamente, fazendo o balanço das correntes nos três 
ramos em paralelo, 


dt | (3.5) 


onde: 


U --= Tensão em V. 
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|, — Corrente do gerador em A, (constante). 


C — Capacidade do condensador em F. 


3.3 Equações de estado do sistema 


Escolhendo para variáveis de estado do sistema 
a temperautra T do termistor e a tensão U do conden- 
sador, O comportamento dinâmico fica descrito pelo sis- 
tema de duas equações diferenciais de primeira ordem 
(3.4) e (3.5) resolvidas em ordem às derivadas, 


| E 
E a = g E [| 
=" O (11,1 + UEM] 


dU 1 
dt Cc to ú 


Os segundos membros destas equações, conside- 
rados como funções das variáveis de estado T e U, são 
fortemente não-lineares. 


(3.6) 


à — AVERIGUAÇÃO DA ESTABILIDADE 


4.1 Linearização das equações 


No estudo das trajectórias definidas no plano de 
fases das variáveis T e U pelo sistema não-linear (3.6), 
a até no cálculo numérico das respectivas solutões, de- 
sempenha um papel fundamental a forma linearizada 
desso mesmo sistema de equações que se obtém consi- 
derando pequenas variações AT e AU das variáveis de 
estado em torno dum ponto fixo de coordenadas T e U 
desse plano. 

O sistema linearizado é da forma, 


AT “a Jo AT 


(4.1) 


LAU ds dA LAG 


onda a matriz de elementos constantes J pq representa 
o jacobiano das funções de T e U que constituem os 
segundos membros das equações (3.6), calculado para 
o ponto considerado. 

De acordo com as equações (3.6), esses elementos 
calculam-so pelas expressões, 


1 ' B 


2 
Jp=— U6 
(4.2) 
1 B 
dg em Ma 
A 
ig tea 
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onde, G é dado por (3.1) e D' se obtém de (3.3) deri- 
vando em ordem a T. 


dD 
P=-—sStPOImTA 


dT (4.3) 


4.2 Pontos de equilibrio 


Igualando a zero os segundos membros das equa- 
ções (3.6), obtém-se o par de equações transcendentes, 


D (T-T,) = U-G(T) 
,=U G(T) (4.4) 


Em relação às trajectórias no plano de fases, estas 
equações definem linhas isóclinas, respectivamente de 
inclinação nula das trajectórias relativamente ao eixo de 
T e do U. Veja-se a Fig. 3. O seu ponto de encontro 
Q define uma situação de equilíbrio, pois que nele se 
anulam simultaneamente as derivadas em ordem ao 
tempo de T e de U, A estabilidade do equilíbrio pode 
ser investigada pela consideração do sistema lineari- 
zado (4.1) referido a esse ponto. 

Qualitativamente, para a situação representada na 
Fig. 3, em que o ponto de equilíbrio se situa no troço 
descendente da isóclina relativa a T estacionário, é ma- 
nifesta a possibilidade desse ponto ser um «foco», se 
so atender aos sentidos de progressão das trajectórias 
(assinalados na figura por setas) e que são determinados 
pelos sinais algébricos assumidos pelos segundos mem- 
bros das equações (3.6) nos quatro sectores definidos 
pelas isóclinas em torno do ponto de equilibrio O. 


LO . e " 
| & ar cê + , 
a a] 
| Ed N ' “e 
f " a ds = Er q 
LED y Mm à Ea 
iadgus Edo ma 
+ f , " 
| 
old] a 
| Ri 
s qm. RE 
j | Rag amd 
| a 
| a 
Db e 


Fig. 3: Linhas Isóclinas no plano de fases 


4.3 Comportamento das trajectórias na vizinhança 
do ponto de equilibrio 


Este estudo faz-se com base na equação caracteris- 
tica da matriz (J) do sistema linearizado (4.1) tomado 


no ponto de equilibrio, 


det Ltd) —- A (1)] =0 
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que é uma equação algébrica do segundo grau em à, 
A2— SpiJJjk + det(J) = O (4.5) 


cujos coeficientes se obtém a partir das expressões (4.2) 
para os elementos de (J), 


[ur 
E 


1 
Sp(J) = 451 + dy — —— (U?6 


B 
(G D' + U2GÊ 7?) 
(4,7) 


Em termos destes coeficientes, as raizes de (4.5) 
são dadas por, 


Sp(J) + V Sp? (J) — 4det(J) | 


(4.8) 


À discussão do comportamento faz-se agora em 
dois passos: 


1.º — Quanto à natureza periódica ou aperiódica 
das soluções de (4.1). 

Esta questão dependa apenas do sinal algébrico 
do discriminante de (4.8). Se se atender à expressão 
(4.7), conclui-se que det(J) é sempre positivo, pelo que 
sa tem: 


a) Sp” (J) < 4det(J). 


As raizes (4.8) formam um par complexo conju- 
gado. Às soluções de (4.1) serão oscilatórias. O ponto 
da equilibrio é um «foco». 


b) Sp” (J) > 4det(J). 


As raizes (4.8) são ambas reais e do mesmo sinal. 
As soluções de (4.1) são aperiódicas. O ponto de equi- 
librio é um anodos». 


2º —. Quanto ao crescimento ou decrescimento das 
soluções com c tempo. 


Esta questão depende unicamente do sinal algé- 
brico de Sp(J). Com efeito, a parte real das raizes (4.8) 
quando são complexas e as próprias raízes quando são 
reais têm o mesmo sinal de Sp(J), pelo que se tem: 


a) Sp(J) > 0. 


Soluções crescentes no tempo. O vonto de equili- 
brio é instável. 


b) SplJ]) = 0, 
Soluções oscilatórias de amnlitude constante. 
c) SpiJ) < 0, 


Soluções decrescentes no tempo. O ponto de equi- 
librio é estável, 
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4.4 Condições de instabilidade 


Para quo a expressão (4.6) seja positiva deve, 
antes de mais, ser positiva a expressão entre parêntesis, 
isto é, 


U?G “73 e ii 


(4.9) 


O significado desta condição pode obter-se consi- 
derando a equação da isóclina a T estacionário (1.º 
equação de 4,4). Diferenciando obtém-se, 


D'- UG — |) dT=2 GU dU 
T2 (4.10) 


A condição (4.9) corresponde a situarmos o ponto de 
funcionamento na região decrescente desta isóclina, 
isto é, onde se tem dU/dT < O como se pode verificar 
em (4.10). Esta é pois a primeira condição necessária 
para se verificar a instabilidade Fig. 3. 

Fixado o ponto de funcionamento, resta somente na 
expressão (4.6) tomar um valor de C suficientemente 
elevado para que a segunda parcela seja inferior à pri- 
meira, ou seja, fazer 

C > Cs (4.11) 
em que C. representa um valor crítico da capacidade 
do condensador. Obtém-se o seu valor igualando (4.6) 
a zero, 


6 k 


Cop =—— = 


o a ' 
(U-G E EEE D') 
E (4.12) 


À condição (4.9) garante um denominador positivo. 


4.5 Periodo de isocronismo das pequenas oscilações 


Quando o ponto de equilíbrio é um foco, o período 
do isocronismo das oscilações pode obter-se a partir 
da parte imaginária das raízes da equação caracteristica. 


Usando (4.8) para Sp” (J) < 4det(J) será, 
2, 2% 
Te —— 
s EImfaki 


1 
VER pd 
V det(J) E Sp“(J) (4:48) 


Para o caso particular de oscilações quase-críticas, 
isto é, para Sp(J) — O tem-se aproximadamente, 


dr 


V det(J) 


Fazendo na expressão (4.7), coerentemente com a 
hipótese, C = C,, dado por (4.12), obtém-se para o 


(4,14) 
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período das oscilações quase-críticas de pequena am- 
plitude, 


(T,) am ( h ) 
rr SS == 
Sum 
onde: 
U“ GB 
à pT7 (4.16) 


Em virtude da condição (4.9) a é sempre maior que 
a unidade nos pontos de funcionamento de interesse. 

O cociento k/D' tem dimensões físicas de tempo 
8 representa a constante de tempo térmica do termistor 
no ponto de funcionamento considerado. O período 
(19) cr é assim proporcional a esta constante de 
tempo. 


o — ESQUEMA EQUIVALENTE INCREMENTAL DO 
TERMISTOR 


Considere-se a equação (3.4) do balanço térmico 
do termistor escrita na forma incremental, na vizinhança 
dum ponto do plano T, U, 


d B 
ka (4T)=(-D' + UÊG 77 ) AT + (206) AU 


(5.1) 
o ainda a relação entre a corrente e a tensão aos termi- 
nais 


|=GU 
também na forma incremental, 

B 
—s )aT + GAU 

T (5.2) 
Introduzindo 4] | definido por, 

at, = 4l — GAU (5.3) 

elimine-se AT na equação (5.1) com o auxílio da (5.2) 


e (5.3). Obtém-se, 


d 
L “o (41, ) + R(4l,) = (AU) 


(5.4) 
cujos coeficientes são dados por, 
B 
T 
L = R = 
2 U26 z 26? 

T* 2 UG 72 

(5.5) 


As equações (5.3) e (5.4) podem imediatamente 
interpretar-se como a descrição do circuito da Fig. 4, 
cujos parâmetros são dados pelas expressões (5.5) em 
termos dos valores de D', U, G e T relativos ao ponto 
de funcionamento considerado. 
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A condição (4.9) que se viu ser necessária para o 
ponto de equilibrio se tornar instável, corresponde, como 
se pode ver em (5.5), a obtermos uma resistência R, 
Fig. 4, com valor negativo. O coeficiente de auto-indução 
L equivalente é sempre positivo. 


8 ád 


Fig. 4: Esquema equivalenta incremental do termistor 
5 — DETERMINAÇÃO EXPERIMENTAL DOS PARAME- 
TROS DO TERMISTOR 
6.1 Medição da temperatura 
Para a verificação experimental da teoria exposta, 
utilizou-se um termistor em forma de vareta ' com as 


dimensões indicadas na Fig. 5 e cujos parâmetros têm 
os seguintes valores nominais: 


mes 
- lá 


Fig. 5: Dimensões do termistor em mm 


Resistência a 25º E R= 150k9 
Constante de temperatura B = 4150 K 
Constante de dissipação a = 13 mwkK” 
Constante de tempo térmica k/S q = 55 s 


fio de cobre 
Bar, // esmaitado 


Fig. 6: Montagem do par termo-eléctrico 


às medidas de temperatura foram feitas usando 
dois pares termo-eléctricos”, um em contacto térmico 
com o termistor e o outro à temperatura ambiente. 

É este conjunto termistor-par, Fig. 6, que é utili- 
zado em todos os ensaios e é a ele que se referem os 
valores da capacidade térmica e das constantes relativas 
à função da dissipação que foram medidos. 


6.2 Caracteristicas eléctricas 


A determinação das constantes características do 
termistor G, e B que figuram em (3.1) foi feita im- 
pondo uma temperatura constante a uniforme ao termis- 
to:, por imersão num banho de óleo, medindo-se a res- 
pectiva condutância em condições de auto-aquecimento 
desprezável. 

No gráfico da Fig. 7 encontram-se representados 
os resultados obtidos para a variação de In(G/G,.| ) 
em função do 1/T, correspondentes à representação em 
gráfico semi-logaritmico da recta, 

G G,, 1 


n— = tn —D- 8 ( — ) 
Gal Gel T (6.1) 


que so cbtém a partir da expressão (3.1). 

Deste resultado experimental conclui-se a validade 
da expressão (3.1) para representar a variação da con- 
dutância com a temperatura no domínio considerado, 
tomando as constantes os valores, 


G, = 4,508 B — 4020 K (6.2) 


o Faliês eaprrmentas 


16 323 2.6 2.9 3.0 31 E 313 24 1 
pod egd + 


Fig. 7: Dependência da condutância relativamente 


à temperatura 


[1h  Termistor da marca Philips n.º 23 226 636 (1,5 W). Veja-se a referência [3] 


[21 Pares Chromel-Alumel em bainha de aço inox isolada. Marca Sodern tipo Thermocoax ZABAeTI, = 43 uV/2C. 
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6.3 Características térmicas 


Para determinar a função de dissipação (3.3), fez- 
-s8 funcionar O termistor no ar em regime de corrente e 
temperatura constantes. No gráfico da Fig. 8 encontra-se 
representada a tensão como função da corrente e no da 
Fig. 9 a potência posta em jogo como função da tem- 
peratura. 

Em regime estacionário teremos, usando a função 
de dissipação na forma (3.3), 


= D=5 (1,1 FE, (6.3) 
Utilizou-se esta expressão para interpolar os pontos 
experimentais, considerando como constantes a determi- 
nar S, OQ e T,, tendo-se obtido: 


S = 10,2mwWwK”! 
OQ = 51,84WKT? (6.4) 
T, = 296,7 
o o Valores experimentais 
1,º235 €C 
[o] 
i E E 3 12 u (má) l 


Fig. B: Característica estática do termistor 


Verifica-se que o valor encontrado para a que 
melhor satisfaz à interpolação é ligeiramente superior, 
cerca de 3,5 “K, ao valor da temperatura do ar ambiente. 

Para determinar a capacidade térmica, realizou-se 
um ensaio em regime dinâmico (dT/dt-*0), utilizando- 
-so a expressão (3.4) na determinação do valor de k, 

Num ensaio de arrefecimento, Fig. 10, determinou- 
-Ss8 a temperatura a cada instante medindo a resistência 
o calculando a temperatura pela expressão (3.1). O valor 
da derivada dT/dt foi calculado numericamente empre- 
gando o método de Lagrange de quarta ordem. Utili- 
zando os valores das constantes da função de dissipação 
já determinados, a substituição de valores na expressão 
(3.4) conduziu a um valor da capacidade térmica: 
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Kk=0773K"! (6.5) 
[5 
[E 
li 
o o Valores experimentais 
12 
TD 
| 
db 
L 
d 
“ 
p* E 
|] 
À | —— ee fee e e 
20 1 39 LO Eq! [7] 70 BO 50 100 (ºC) 
- i = —— 
150 de nb 30 330 350 3150 350 IO PK) T 


Fig. 9: Dissipação em função da temperatura 


6.4 Valores relativos a um ponto de funcionamento 
típico 


Para o estudo do comportamento dinâmico do sis- 
tema condensador termistor foi escolhido o ponto de 
funcionamento da característica estática correspondente 
al=4maA, Fig. 8. Obtém-se nesse ponto, 


U=129V T = 338 K (6.6) 
Os valores dos parâmetros do esquema equivalente 

da Fig. 4 podem agora ser obtidos para este ponto de 

funcionamento, fazendo uso das expressões (3.1) e (5.5) 

a são, 

G =3145 L = 685 kH 


R = — 3,380 kQ (6.7) 


Bo 


Tn 20 PC 


0 20 £0 Eq 80 00 Ro 10 


150 (s)t 
Fig. 10: Curva de arrefecimento 


7. — SIMULAÇÃO NUMÉRICA 


Utilizou-se um método preditor-corrector para a 
integração numérica das equações diferenciais do sis- 
tema. O método baseia-se no emprego, na fase de cor- 
recção, de fórmulas lineares multipasso, concretamente 
das fórmulas de Gear [4], [5] de ordem um a quatro. 
A informação quanto ao passado da solução é retida sob 
a forma de coeficientes de polinómios intepoladores, um 
para cada variável da solução, de grau q igual à ordem 
da fórmula correctora empregada e que passam pelos 
q + ? últimos pontos calculados, A predição faz-se por 
simples extrapolação para o novo ponto da solução a 
calcular, usando os referidos polinómios. 

O método preditor-corrector é utilizado como parte 
integrante dum algoritmo adaptativo, (algoritmo de Gear 
[5]), que estabelece automaticamente o passo h e a or- 
dem q da fórmula a empregar, tendo em vista maximizar 
o passo h mantendo todavia os erros locais abaixo dum 
certo limite especificado pelo utilizador. O arranque a 
partir das condições iniciais do problema processa-se 
automaticamente pelo emprego da fórmula de primeira 
ordem, 

Na aplicação do algoritmo, as equações pertinentes 
do problema intervém apenas na fase de correcção. Re- 
presentando o sistema de equações a integrar abreviada- 
mente na forma, 


(7.1) 


onde y é o vector das variáveis de estado, T e U neste 
caso, a F, é o vector das funções dessas variáveis que 
constituem os segundos membros das equações (3.6). 
A operação de correcção consiste em determinar incre- 
mentos Ay e Ay a adicionar aos valores de y e w 
obtidos na predição, tais que as equações (7.1) sejam 
satisfeitas no ponto considerado sob a condição de ser, 


Ay =boh dy (7.2) 
onde h é o incremento da variável tempo e Do é uma 
constante dependente da ordem qg da fórmula utilizada 
[4]. 

Quando as funções F y são não-lineares, a determi- 
nação dos incrementos consegue-se por um processo 
iterativo segundo o método de Newton-Raphson que se 
baseia na forma linearizada do sistema (7.1), 


vo F, + ay —(J),ày=0 


onda (JJ), representa a matriz do jacobiano das funções 
F, cujos elementos são dados pelas expressões (4.2) e 
(1) representa a matriz unidade. Conjugando esta expres- 
são com (7.2) obtém-se um sistema de equações linea- 
res para os incrementos Ay, 


— (1)-— ayv=—|y' — F 
o li, | ar ai 


O cálculo das correcções inicia-ss tomando para 
y e y os valores obtidos na predição, com os quais se 


(1) Marca Tektronix, modelo TEK-31. 


calculam as funções F, e (JJ, que permitem estabe- 
lecet o sistema (7.3). Obtidos os incrementos Ay por 
solução de (7.3), inicia-se nova iteração partindo de va- 
lores actualizados de y e y' resultantes de fazer 


gr Ay 
e atendendo a (7.2) 
yo+ ay(t/b h)>y 


O processo repete-se até à convergência, isto é, 
até que, segundo o critério de tolerância de erros adop- 
tado, os incrementos Ay se tornem desprezáveis. Nor- 
malmento são suficientes duas iterações. 

O algoritmo foi programado em termos gerais para 
uso numa máquina de calcular de secretária”, requerendo 
apenas para cada problema específico a programação 
das fórmulas para o cálculo de F, e (J),, do problema. 
Na integração das equações (3.6) cada passo de inte- 
gração requer em média cerca de 155. 


8 — ESTUDO EXPERIMENTAL DO SISTEMA 
CONDENSADOR-TERMISTOR 


Não se dispondo de condensadores com a capa- 
cidado necessária (= 10' «uF) para os valores da tensão 
exigidos (= 200 V), simulou-se com o auxílio da mon- 
tagem da Fig. 11, empregando amplificadores operacio- 
nais e um triodo (afim de acomodar os elevados valores 
da tensão), não só o efeito do condensador como o 
próprio gerador ideal de corrente do esquema da Fig. 2. 

Com efeito, em relação ao par de terminais AB da 
Fig 11, onde fica ligado o termistor, a montagem define 
uma relação corrente-tensão que, admitindo os amplifi- 
cadores como ideais, é da forma, 


dU 
=), = 


he á 


(8.1) 


idêntica à equação (3.5) descritiva do circuito eléctrico 
da Fig. 2, se atendermos a que i representa o valor ins- 
tantâneo da corrente no termistor. 

Os coeficientes dae (8.1) podem definir-se a partir 
dos elementos que figuram no esquema da Fig. 11 pelas 
axpressões, 


(8.2) 


(8.3) 


Na realização prática, o factor multiplicativo da 
capacidade C, na expressão (8.2) era de cerca de 170, 
tendo o valor da capacidade C simulada sido medido 
num ensaio de «descarga» sobre uma resistência cons- 
tanta de valor conhecido ligada entre A e B, Fig. 11 
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Re=5kQ Ru=z 2 MQ 
R= 1 MQ Ru= 1 M2 
R:a= 1 Mº R = 10 k02 
R,=.5 MQ 
Fig. 11: Montagem para simulação do condensador e do 


gerador ideal de corrente 


A determinação da temperatura do termistor foi 
feita por intermédio do valor da sua condutância, como 
se descreveu em 6.3. Em regime dinâmico, a tempera- 
tura obtida através do par termo-eléctrico não corres- 
ponde ao valor instantâneo da temperautra do termistor, 
em consequência da resistência térmica do contacto 
par-termistor e inércia têrmica do par. Praticamente 
observa-se um atraso no tempo da indicação do par re- 
lativamente à temperatura instantânea. Ver Fig. 12. 

A verificação experimental do comportamento di- 
nâmico do sistema condensador-termistor foi feita fi- 
xando |, = 4mAÃ e C = 9880 uF. Utilizando os valores 
anteriormente calculados (6.4), (6.5) e (6.6) na expres- 
são (4.12) obtém-se para este ponto de funcionamento 

U 
(V) 


140 


120 


100. 


801 | 


60 | 


20 30 LO 50 60 70 


uma capacidade crítica C,., = 6490 4F pelo que se 
trata duma situação instável C = 1,52 C,, . No domínio 
das pequenas oscilações, pode prever-se um periodo 
(4.13) de 547 s e um incremento da amplitude 1,58 
vezes por período. 
PEL mr) 

1.8 
120 —— Temperatura (res) 


LE — — Temperatura (par) 


-=-— Potencia 


109 200 300 100 S00 500 700 809 (5) 4 


Medida dinâmica da temperatura. Comparação 
da resultados 


Fig. 12: 


Na Fig. 13 estão traçadas as trajectórias no plano 
de fases calculadas para as condições das experiências 
utilizando o método de simulação numérica descrito em 
7. o sobre elas os pontos que foram obtidos experimen- 


Valores calculados 


— — — 


* x Valores médidos 


- 8.0] 
EE O = 4 mA 


Bo 90 100 no 0 30 (ºC) 1 


Fig. 13 — Comparação das trajectórias calculadas com os valores experimentados 


— Walcrts caltujados 


=. .Valgret meddos 


ii o 
500 Cafe 00 1200 DO 1660 (ue) + 


f a 


e Dr er ste (— e e e e e e 


Fo Lua ELU 200 1000 1200 1400 1600 (9) + 


Fig. 14: Comparação dos valores calculados e experimentais das 
variáveis de estado como função do tempo 


talmente. Verifica-se efectivamente a existência dum 
ciclo limita com a forma e dimensões previstas, conver- 
gindo as trajectórias para ela quer as condições iniciais 
sa situem no interior quer no exterior do ciclo. 

Na Fig. 14 apresentam-se os mesmos resultados 
teóricos e experimentais na sua dependência do tempo 
a no que respeita ao regime periódico final, (ciclo li- 
mita da Fig. 13). 4 concordância das formas de curva 
calculadas com os pontos experimentais é muito boa. 


9 — RESULTADOS OBTIDOS POR SIMULAÇÃO 
NUMERICA 


Para além da situação investigada em 8., utilizou-se 
a simulação numérica para averiguar o comportamento 
do sistema para outros valores da capacidade, mantendo 
as restantes condições de funcionamento. Consideraram- 
-so dois casos extremos, C = 0,5 €.., eC=2C,,. 

O primeiro caso corresponde a uma situação está- 
vel. Largando o sistema com o termistor à temperatura 
ambiente e o condensador descarregado, ele vem a atin- 
gir a situação de equilíbrio mas por intermédio de um 
processo oscilatório amortecido, Fig. 15. Durante o pro- 
cesso de evolução para o equilibrio o sistema passa por 
valores da temperatura e da tensão, particularmente 
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desta última, consideravelmente superiores aos de equi- 
librio. 

No segundo caso o ponto de equilibrio é instável 
o as trajectórias no plano de fases tendem para um 
ciclo limite, Fig. 16. Dada a situação fortemente sobre- 
-crítica o ciclo é muito mais amplo que o da Fig. 13, 
particularmente segundo o eixo da temperatura que 
atingo valores que seriam destrutivos para o termistor 
usado nas experiências. 

Na Fig. 17 comparam-se os dois casos sob o ponto 
do vista da evolução das variáveis de estado com o 
tempo e para idênticas condições iniciais. Para além 
do carácter distinto do regime final, constante num 
caso, periódico no outro, ressalta a diferença dos perio- 
dos de oscilação quo se encontram sensivelmente na 
mesma relação quo as raízes quadradas dos valores das 
capacidades, cerca de 1 para 2 neste caso. 

Em relação aos valores de equilíbrio nota-se que 
a sobre-elevação da tensão predomina no caso da capa- 
cidade menor e a sobre-elevação da temperatura no 
caso da capacidade mais elevada. Neste último caso a 
curva da tensão aproxima-se de forma dum dente de 
serra como é típico de oscilações de relaxação, coinci- 
dindo a queda rápida da tensão com o pico de tempera- 
tura do termistor. 

Os gráficos da Fig. 17 permitem uma interpretação 
física simples do processo oscilatório, neste caso ex- 
tremo das capacidades elevadas: 

Durante a fase crescente da tensão o termistor está 
frio, a sua resistência é elevada, assiste-se assim prati- 
camento à carga do condensador pela cerrente cons- 
tanta do gerador. Dal o crescimento da tensão ser pra- 
ticamento linear com o tempo. Pode constatar-se tam- 
bém, especialmente na fase de arranque, que a veloci- 
dade de crescimento da tensão é inversa do valor da 
capacidade. ' 

Enquanto o condensador carrega, o progressivo au- 
mento da dissipação no termistor eleva-lhe a tempera- 
tura baixando-lha a resistência até ao ponto em que, 
para a tensão existente no condensador, toda a corrente 
do gerador circula agora pelo termistor, parando assim 
a tensão de crescer. Daqui em diante o progressivo au- 
mento da temperatura do termistor tende a baixar a 
tensão e o condensador começa a descarregar-se sobre 
o termistor. Inicia-se assim um processo comulativo em 
avalanche, pois que agora além da corrente do gerador, 
contribui para o aquecimento do termistor a corrente de 
descarga do condensador produzindo-se um rápido au- 
mento da temperatura e rápida descarga do condensaodr 
através da baixa resistência do termistor. À diminuição 
da tensão a consequente diminuição da dissipação no 
termistor permite o seu arrefecimento e o início de um 
novo ciclo, 

Finalmente na Fig. 18 mostra-se o nascimento es- 
pontâneo das oscilações quando o sistema é largado na 
proximidade do ponto de equilibrio no caso instável e 
qua corresponde à espiral interior do diagrama da Fig. 
16. As oscilações de pequena amplitude são pratica- 
mente sinusoidais, afirmando-se progressivamente o 
carácter fortemente não-linear à medida que a amplitude 
cresce. Notar que a estabilização do regime periódico 
exigiria tempos superiores a uma hora. 


TÉCNICA 431 


(V) 
200- | 
1 
7 C=0.5C 
, cr 
180 ! 
| N 
| N 
1607 À 
: 
1 
1404 > 
e qu e um É N 
= E me 
+ ' O ma 
120 Fá SE 
, 4 Ea dL.o 
/ 4 ai O dt 
Fr a 
x E A 
100 a N = 
" 
" 
) ; " 
4 
Í " 
BO E 
| | q 
| a 
| | E 
50 | de 
| e 
| E 
Í AS 
| Í cia 
0 | cá ea 
“ Mn 
dt lg=b má 
| 
20 
| 
20 30 LO 50 50 70 Bo 90 100 (ºC) T 
Fig. 15: Trajectória no plano de fases. Caso sub-critico 
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Fig. 17: Evolução das variáveis de estado no tempo. Comparação dos casos sub e sobre-crítico 
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Fig. 18: Evolução das variáveis de estado no tempo. Arranque das oscilações 


10 — CONCLUSÕES 


O sistema condensador-termistor é sede de oscila- 
ções auto-excitadas quando o ponto de funcionamento 
em repouso determinado pelo gerador ideal de corrente 
que alimenta o sistema se situa no troço descendente da 
característica estática corrente-tensão do termistor, e a 
capacidado do condensador excede um certo valor cri- 
tio. À não-linearidade do termistor assegura a existência 
dum ciclo limite estável. O período das oscilações é 
essencialmente condicionado pelo valor da capacidade 
o pela inércia térmica do termistor. 

Verifica-se que as propriedades do termistor podem 
ser descritas adequadamente por uma lei exponencial 
para a dependência da condutância em relação à tem- 
peratura e por uma lei de dissipação função da diferença 
de temperatura entre o termistor e o ar ambiente 
envolvendo um termo linear e um termo quadrá- 
tico de correcção. Esta descrição conduz, no domínio 
incremental, a um esquema eléctrico equivalente de tipo 
indutivo onde surge uma resistência que toma valores 
negativos para pontos de funcionamento no troço decres- 
cente da caracteristica estática. 

O comportamento dinâmico do sistema condensa- 
dor-termistor calculado através da simulação numérica 
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do modelo teórico proposto é confirmado com bastante 
rigor pelas medidas efectuadas numa montagem expe- 
rimental., 
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RESUMO 


Recordam-se as duas revoluções sofridas pela En- 
genharia de Estruturas, provocadas respectivamente pela 
utilização dos princípios da Mecânica ,no século XVIII, 
e dos computadores digitais em meados do século XX, 
O desenvolvimento de técnicas de discretização, como a 
dos elementos finitos, veio permitir a resolução de nu- 
merosos problemas. Sente-se no entanto a necessidade 
de uma profunda modificação dos métodos tradicionais 
que torne possivel uma mais eficaz exploração das po- 
tencialidades das máquinas. 


1 — PARA ALÉM DA GEOMETRIA 
1.1 — Instrumentos matemáticos anteriores ao séc. XVIII 


Uma das interrogações que surgem no espírito ao 
contemplar certas grandes estruturas, antigas, medie- 
vais, renascentistas ou barrocas, é a de quais pode- 
rão ter sido os métodos empregados no projecto dessas 
cobras, já que os instrumentos científicos básicos de que 
se serve o moderno projectista são de utilização relati- 
vamente recente, 

Sabe-se que se usavam regras empíricas de natureza 
geométrica [1] no projecto das fundações, pilares, abó- 
badas, etc., e que essas regras, que procuravam garantir 
a estabilidada das construções, se empregavam junta- 
mente com outras de significado meramente estético. 

Podem servir de exemplo as regras de Alberti para 
pontes de pedra, as quais são citadas num livro [2] 
impresso em Florença em 1485: «A largura dos pilares 
devo ser 1/4 da altura da ponte; o vão livre do arco 
não deve ser maior que 6 vezes nem menor que 4 vezes 
a largura dos pilares; a espessura das aduelas não deve 
ser menor que 1/10 do vão». 

Os construtores seguiam tais princípios da arte 
um pouco como hoje os engenheiros seguem os regula- 
mentos, À grande diferença é que estes constituem o 
produto final de todo um conjunto de investigações, en- 
quanto que às regras antigas faltava a base científica. 

Considerações de ordem geométrica estiveram 
ainda presentes, como é evidente, na brilhante resolu- 
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SUMMARY 


The two revolutions suffered by Structural Engine- 
ering, respectively in the XVIII th century, due to the 
introduction of the principles of Mechanics, and in the 
XX th century, due to the advent of computers, are first 
remembered. Attention is then called to the fact that a 
fully efficient exploration of the capacity of computers 
can be possible only if a deeper change of traditional 
methods is achieved. 


ção bizantina do problema da cúpula esférica sobre 
planta quadrada e na não menos brilhante resolução 
gótica do da cobertura de um espaço rectangular por 
uma abóbada apoiada exclusivamente nos quatro cantos, 
recorrendo à chamada cruzeta de ogivas [3,4]. 

Acreditava-se que uma solução geometricamente 
eleganto era uma solução estruturalmente correcta. 

Tem-se hoje bem consciência de que, embora uma 
boa solução estrutural tenha de ser geometricamente ele- 
gante, nem todas as soluções geometricamente elegantes 
são boas soluções estruturais. Não é por exemplo a geo- 
metria que nos ensina que a forma parabólica é, para um 
arco de betão, mais lógica que a forma circular, embora 
esta o possa ser para o arco de pedra porque permite 
que as aduelas sejam todas iguais (raciocínio puramente 
geométrico). À geometria também pouco pode dizer-nos 
sobre as dimensões a dar aos elementos de uma estru- 
tura, e menos ainda sobre o modo como essas dimen- 
sões devem variar com o material utilizado 


1.2 — O nascimento da moderna Engenharia 
de Estruturas 


à engenharia moderna nasceu quando os princípios 
da Mecânica e da Física passaram a ser aplicados na 
resolução dos seus problemas, 

A primeira tentativa de aplicação da Mecânica na 
construção fez-se por volta de 1742, quando, por ordem 
de Bento XIV, os matemáticos Le Seur, Jacquier e Bos- 


(*) Conferência preparada para as adiadas 4.º jornadas Luso-Brasileiras de Engenharia Civil. 
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covich tentaram a análise estrutural da cúpula de 5, Pe- 
dro do Roma com o fim de estudarem um processo para 
a eliminação das fendas que nela tinham aparecido [1]. 

O referido estudo levou à introdução de cinco anéis 
metálicos que, cintando a base da cúpula, absorveram 
os respectivos impulsos. Foram postas de lado as de- 
sastrosas soluções apresentadas noutros pareceres, as 
quais envolviam, umas, o enchimento das caixas de es- 
cada, outras a remoção da lanterna, outras ainda a 
construção de quatro contrafortes colossais. 

O parecer de Le Seur, Jacquier e Boscovich foi 
muito atacado. 

Disse-se por exemplo que os cálculos não podiam 
estar certos porque, se estivessem, já a cúpula deveria 
ter caido. E esta critica era justa porque na realidade os 
cálculos estavam errados. 

Muito menos razão tinham aqueles que os atacaram 
com base na ilegitimidade da aplicação da Mecânica 
e da Fisica: ter feito essa aplicação foi justamente o 
grande mérito dos três matemáticos. 

A utilização dos princípios da Mecânica determinou 
o aparecimento da Análise Estrutural e desencadeou 
uma revolução que transformou totalmente as condições 
em que os construtores exerciam a sua actividade. Data 
dai de resto a dramática separação entre engenheiros e 
arquitectos. 

Esta revolução, a que chamaremos a primeira, teve 
certamente inconvenientes. 

Passou-se por um lado a só construir aquilo que se 
sabia calcular. Perdeu-so por outro lado intuição e ima- 
ginação. E o facto de os construtores se terem dividido 
em arquitectos e engenheiros teve nefastas implicações 
na arquitectura do século XIX. 

As vantagens foram porém muito grandes. 

O aparecimento da Análise Estrutural permitiu com 
efeito a obtenção de soluções mais económicas, a demo- 
cratização da técnica e a facilidade de extrapolação com 
base em teorias científicas. E só assim a Engenharia 
de Estruturas pôde responder aos grandes desafios que 
lha foram postos desde os princípios da Idade Contem- 
porânea, desafios esses que resultaram da necessidade 
da construir muito com recursos limitados, de construir 
depressa e com segurança, de se adaptar constantemente 
a novas situações e, finalmente, de aumentar largamente 
o número de projectistas, cometendo a homens vulgares 


8 não só a génios o encargo de projectar as grandes 
construções, 


1.3 — O impacto dos computadores 


Originada pelo aparecimento dos computadores, 
deu-se na Engenharia de Estruturas, por volta dos mea- 
dos do século XX, aquilo a que pode chamar-se a se- 
gunda revolução, 

Esta segunda revolução foi em grande parte deter- 
minada pelas novas e extraordinárias possibilidades de 
cálculo, as quais, permitindo que a Técnica se servisse 
mais largamente dos instrumentos científicos existentes, 
levaram os engenheiros a penetrar em domínios que 
antes consideravam fora dos seus interesses. 
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Assim, por exemplo, antes do aparecimento dos 
computadores, a determinação de tensões, deformações 
e deslocamentos em corpos elásticos solicitados por 
acções exteriores era objecto de duas ciências distintas: 
a Teoria da Elasticidade, caracterizada pelo rigor e cul- 
tivada pelos matemáticos, e a Resistência de Materiais, 
subordinada à eficácia e desenvolvida pelos engenheiros. 

Ora, do impacto dos computadores, que vieram 
permitir a resolução numérica das equações diferenciais 
da Teoria da Elasticidade, e, por isso, a sua utilização 
corrente pelos engenheiros, resultou a sintese de ambas 
numa única ciência, a Mecânica dos Sólidos, pela qual 
so interessam simultaneamente hope em dia os enge- 
nheiros e os matemáticos. 

A segunda revolução teve assim a virtude de eli- 
minar o inconveniente originado pela primeira de se 
ter passado a construir exclusivamente aquilo que se 
sabia calcular. Tal inconveniente deixou efectivamente de 
o ser desde que as equações diferenciais dos vários mo- 
delos estruturais puderam ser objecto de tratamento nu- 
mérico, isto é, desde que todas as estruturas puderam 
ser calculadas. 

Aos métodos de discretização que permitiram tão 
extraordinária conquista será dedicada a segunda parte 
deste artigo. Observamos contudo qua para essa con- 
quista muito contribuiu também o desenvolvimento dos 
métodos experimentais, aos quais só não fazemos refe- 
rência mais desenvolvida porque não desejamos afastar- 
-nos da nossa especialidade. 

Os engenheiros não se contentaram porém com a 
utilização dos conhecimentos científicos existentes, já 
que sob a sua solicitação se empreenderam e empre- 
endem numerosas investigações destinadas ao estabele- 
cimento de novas teorias capazes de simularem mais per- 
feitamente o mundo real e de resolverem problemas mais 
próximos das necessidades da Técnica. Alguma coisa se 
clirá quanto a este ponto na terceira parte, 

Finalmente, as consequências do aparecimento dos 
computadores tendem a ser ainda mais profundas, dado 
que a Técnica, não se limitando a buscar na Ciência ins- 
trumentos para resolver problemas mais ou menos clás- 
sicos, tende a reformular a sua própria problemática [5]. 

E o que parece extremamente animador é que tal 
reformulação é feita, como se verá na última narte deste 
trabalho, no sentido de libertar o engenheiro do segundo 
a grande inconveniente da primeira revolução: o da 
perda de sensibilidade e de imaginação. 


2 -— PARA ALÉM DA ESTÁTICA 
21 — O problema da discretização 


Antes do aparecimento dos computadores digitais, 
a Teoria das Estruturas incluia unicamente o estudo das 
estruturas reticuladas, estudo esse a que a necessidade 
de evitar a resolução de sistemas de muitas equações 
a muitas incógnitas emprestava dificuldades e deselegân- 
cias que os computadores, mais ou menos indirecta- 
mente, muito ajudaram a fazer desaparecer. 

A partir da segunda metade da década de quarenta, 
a complexidade progressiva dos cálculos (sobretudo na 
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Aeronáutica), as perspectivas revolucionárias que se 
abriram à resolução de grandes sistemas de equações e 
a necessidade de formulações que facilitassem a automa- 
tização levaram à introdução do formalismo matricial na 
Teoria das Estruturas [6]. 


Dez anos mais tarde, isto é, nos anos cinquenta, O 
problema punha-se já em termos de estender essas for- 
mulações à análise de estruturas não reticuladas, 


A dificuldade estava toda em que somente as es- 
truturas reticuladas são naturalmente discretas, isto é, 
só elas gozam do privilégio de, enquanto decomponiveis 
em elementos idealizáveis como peças lineares, O seu 
estado de deformação poder ser descrito exactamente 
em termos de um número finito de parâmetros e, por- 
tanto, de a sua análise ser redutível, de modo exacto, à 
resolução de sistemas de equações algébricas [7]. 


Já o mesmo não acontece efectivamente com as 
estruturas não reticuladas, que são regidas por equações 
diferenciais com derivadas parciais e cuja discretização 
só pode ser conseguida mediante a introdução de hipá- 
teses aproximadas. 


É claro que o problema da resolução numérica de 
equações diferenciais tinha sido abordado muito antes 
da década de cinquenta, quer através de métodos de 
diferenças finitas, quer exprimindo a solução em termos 
do um número finito de parâmetros, em regra coeficien- 
tes de uma combinação linear de funções, e determi- 
nando esses parâmetros por diversas técnicas dos quais 
a de Ritz é talvez a mais notável, quer formulando o 
problema em termos de sistemas de equações integrais 
e promovendo a respectiva resolução numérica, método 
esto quo mais tarde daria origem à classe mais geral 
dos métodos integrais. 

No ano de 1956 foi ho entanto publicado um ar- 
tigo [8] de importância singular na História da Análise 
Estrutural, no qual se descrevia pela primeira vez a apli- 
cação de um método, denominado dos elementos finitos, 
à resolução de problemas de elasticidade plana. 

O método não era novo, a bem dizer, dado que já 
em 1943 Courant tinha publicado um trabalho [9] sobre 
a resolução de problemas de torção utilizando uma téc- 
nica que, no fundo, era a dos elementos finitos. Também 
Synge [10] trabalhara na mesma técnica dede 1952, 
a só o facto de os tempos não estarem maduros, pela 
deficiência dos meios de cálculo então utilizados, ex- 
plica que a obra dos dois notáveis matemáticos não 
tenha dado origem imediata à explosão que alguns anos 
mais tarde se verificaria, 

Quais as vantagens da nova técnica? 

Tal como a das diferenças finitas, ela pertencia 
à classe dos métodos, sempre sedutores para o enge- 
nheiro, que, embora com sacrifício da precisão, apresen- 
tam uma extrema generalidade, 


Oferecia no entanto, relativamenta à técnica das 
diferenças finitas, outras vantagens que iriam permitir- 
“lho destronar esta última. 

Em primeiro lugar, libertava o utilizador da tirania 


do uma rede regular. Permita, em segundo lugar, mais 
fácil tratamento das condições de fronteira, 

Estas duas vantagens que, como veremos, eram mais 
aparentes que reais, não teriam sido no entanto tão im- 
portantes como uma terceira para a popularidade que 
o método dos elementos finitos veio a alcançar junto 
dos engenheiros. E essa foi a de ter permitido uma gene- 
ralização fácil dos conceitos e resultados da teoria das 
estruturas reticuladas, 

Os engenheiros encontraram portanto um método 
próximo daqueles com que estavam habituados a traba- 
lhar. Mas não só, já que a tal vantagem, de ordem 
principalmente psicológica, sa somava a facilidade da 
automatização, aproveitando todo o trabalho já desenvol- 
vido no sentido da formulação matricial da Análise 
Estrutural, 

Por outro lado, no plano mais teórico da Mecânica 
dos Sólidos, a técnica dos elementos finitos trazia a pos- 
sibilidade de construir um modelo discreto inteiramente 
análogo aos modelos contínuos tri, bi e unidimensio- 
nais, modelo discreto esse que, chamando a atenção 
para o papel extremamente importante desempenhado 
pelo conceito de convergência [11], iria lançar nova luz 
sobre questões como a dos fundamentos da teoria das 
cascas [12], permitindo a elaboração de uma teoria 
matemática das estruturas e, talvez mesmo, a de uma 
mecânica dos meios discretos. 


2.2 — O método dos elementos finitos como técnica de 
interpolação * 


Considere-se o problema clássico da Análise Nu- 
mérica que consiste em determinar uma função f que 


tome em abcissas dadas, x, , valores fixados f, , isto 
é, que verifique as condições: 
fix )=f,(i=1,... n) (1) 


Diz-se que f (x) interpola os valores fixados para a 
função nas abcissas x;,8 o problema correspondente é o 
da interpolação a uma dimensão no intervalo (x, x,). 

Dois tipos de soluções têm sido preconizados para 
a resolução deste problema, 

No primeiro tipo, clássico, a função f (x) tem expres- 
são analítica, em geral polinomial, válida sobre todo o in- 


tervalo (x, X,). € pertence à classe Cx axa). isto é tanto 
a função como todas as suas derivadas são contínuas no 
intervalo, 

No segundo tipo, a função f (x) é definida seccional- 
mente, isto é, tem expressões analíticas distintas em di- 
ferentes intervalos (x, , x, +, ), sendo assegurada nas 
extremidades dos intervalos a sua continuidade e a das 
suas derivadas até à ordem r. Tem-se então o que se 
chama uma função «spline» (**) da classe C'(x,, x, ). 

Quando se comparam vários processos numéricos 
aplicados com um mínimo de competência, observa-se en 
geral a regra de ganhar por um lado o que se perde por 


(*) Relativamente a esta secção e à seguinte, chama-se a atenção para [13]. 


(**) «Spline» é uma palavra inglesa que designa uma régua flexivel de madeira ou metal utilizada para desenhar curvas. 
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outro. Tudo está portanto em adaptar o processo ao 
problema a resolver, perdendo o que não faz falta para 
ganhar o que pode ter utilidade. 

No caso presente, a vantagem das soluções do se- 
gundo tipo é evidente, já que permitem uma muito maior 
liberdade na construção das funções, O que com elas 
se perde não é por outro lado muito significativo porque 
não se tem normalmente necessidade de que sejam con- 
tinuas as derivadas além de certa ordem. Dal a grande 
superioridade destes métodos. 

A duas dimensões podemos ter também soluções 
dos dois tipos para o problema que consiste em cons- 
truir uma função que tome valores fixados em pontos, 
a que daremos o nome de nós, de um domínio bidimen- 
sional. Corresponde naturalmente o segundo tino à sub- 
divisão do domínio em subdomínios poligonais contendo 
os nós nas suas fronteiras, em regra como vértices, e à 
utilização de funções interpoladoras com expressões ana- 
líticas distintas em cada um desses subdomínios, 

A construção dos «splines» pode tornar-se no en- 
tanto extremamente dificil, sobretudo se os dominios fo- 
rem triangulares e não rectangulares, e se se quiser 
impor a continuidade não só da função como de deriva- 
das de qualquer ordem, inclusivamente da primeira. 

Pode utilizar-se no entanto a técnica dos elementos 
finitos que se distingue da dos «splines» porque a con- 
tinuidade é imposta nos nós e não necessariamente em 
todos os pontos das fronteiras entre elementos. 

Observa-se que no caso unidimensional não faz sen- 
tido distinguir entre as duas técnicas, porque a fronteira 
se reduz a um único ponto, que é assim simultaneamente, 
um nó, À distinção surge precisamente a partir das duas 
dimensões, 

Um elemento finito é então um subdomínio fechado 
associado a um espaço linear de funções definidas sobre 
esse subdomínio, 

Tal espaço linear pode definir-se dando a função 
f [x) como função linear dos valores interpolados em 
cada elemento, isto é, fazendo 


xe AS (2) 


o 


f (x) - qe T (x) fe para 


em que A“ designa um subdomínio genérico e o índice 
e se refere a tal subdominio. f . representa o vector dos 
valores que a função toma nos nós do elemento, isto é, 
dos valores interpolados no elemento e. 

As funções q; (x), elementos do vector q“ 
tomam valores nulos em todos os nós do elemento 
excepto no nó i. Dizem-se funções de forma («shape 
functions») [14] associadas com o elemento e. 

Denominando nó do domínio A qualquer dos pontos 
do domínio que seja nó de um ou mais subdomínios Ae, 
podemos dar à função f (x) uma expressão válida para 
todo o domínio A, desde que se considere o vector f dos 
valores interpolados em todo o domínio. Tem-se então: 


fh)=eTt=o,f (3) 


em que q, é uma função de forma associada com o 
dominio e não já com um elemento, 

A três dimensões, o problema põe-se praticamente 
como a duas dimensões, a única diferença residindo em 
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que, sendo muito mais elevado o número de valores a 
interpolar, valerá a pena utilizar elementos mais comple- 
xos, isto é, abrangendo maior número de nós, ou consti- 
tuidos à custa de elementos mais simples. 

O emprego de coordenadas naturais [14] dos ele- 
mentos, além de oferecer outras vantagens, pode facilitar 
muito a determinação directa de funções de forma. En- 
tende-se por sistema de coordenadas naturais um sistema 
de coordenadas locais (isto é, relativo a um determi- 
nado elemento e não ao sistema de elementos) que per- 
mite especificar a posição de um ponto dentro do ele- 
mento por um conjunto de grandezas adimensionais 
cujos valores variam entre zero e um, valores extremos 
estes que geralmente se atingem nos nós. 

As chamadas coordenadas de área ou volume são 
as coordenadas naturais utilizadas respectivamente para 
os triângulos e tetraedros. Para os rectângulos, as coor- 
denadas naturais são as cartesianas com origem no bari- 
centro e eixos paralelos aos lados, 

Verificando-so que há dificuldada em cobrir con- 
venientemente domínios de fronteira curva com elemen- 
tos da lados rectilineos sem exagerar o número destes, 
há hoja em dia tendência para utilizar elementos com 
lados curvos [15], os quais podem ser vistos como 
resultantes da deformação de subdomínios com lados 
rectos por uma transformação de coordenadas do tipo 


x=pel (E) xe; 


ams 


y=pelt)ye;z=qel (E) ze 
(4,5,6) 
em que qt representa um vector de funções de forma. 

Se também para a função f se tiver a expressão 
traduzida pela equação (2), o elemento diz-se isopara- 
métrico. Será sobreparamétrico se forem necessários 
mais nós para definir a geometria que a função, sub- 
paramétrico no caso contrário, que é o caso geral dos 
elementos de lados rectilineos. 

Os elementos isoparamétricos, cuja popularidade 
tem vindo a acentuar-se, gozam de propriedades extre- 
mamente interessantes. 

Assim, se as funções de forma forem tais que se 
anulem sobra os lados (ou faces, no caso tridimensio- 
nal) que não contenham o nó a que se referem, a con- 
tinuidade nos nós implica a continuidade total. Obtemos, 
por outras palavras, subdomínios que se conformam uns 
com os outros, sem rasgos nem sobreposições, e funções 
interpoladoras da classe Cº. 

Uma segunda propriedade, que é fundamental para 
garantir a convergência, exprime-se por um teorema que 
afirma que, se as funções de forma forem tais que 
e ve = 1,a função poderá ter dentro de cada elemento 
uma variação linear arbitrária. 


2.3-=- 0 método dos elementos finitos como método 
de aproximação 


Seja de novo um domínio A com fronteira B, e 
considere-se uma equação diferencial 


Lf=g em À (7) 
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com condições de fronteira 
Cf=h em B (8) 


A resolução exacta de tal equação será problema 
muito difícil. Pode ser no entanto muito mais fácil de- 
terminar uma solução suficientemente próxima da exacta 
num espaço E. com número finito, D, de dimensões. 

A técnica dos elementos finitos intervém justamente 
na construção do espaço E, - Supor-se-á Dois no que 
se segue que os elementos de E,, são as funções f(x) 
definidas pela equação (3). 

Várias técnicas podem ser utilizadas para a deter- 
minação da solução aproximada [13]. 

A técnica mais óbvia é a das colocações, que con- 
sisto em seleccionar a função do espaço Ep que satis- 
faça a equação (7) ou as condições de fronteira (8) 
num número D de vontos respectivamente do domínio 
ou da sua fronteira. 

É-sa assim conduzido ao sistema de equações 


(E qmlfm=9, para P, cá (9) 


(C ymlntm=— h, para P, «B (10) 


em que a simboliza o nó de ordem n. 
Outro método possível é o dos minimos quadrados, 
que corresponde a procurar f = Ep tal que o integral 


=D [ut-graat+ [tct-ni dp -0 
 AÉ B 


(11) 


seja minimo. 
Introduzindo [3] tem-se 


= D) | [Lp,)ty— 9] dat + 


fre pa) fa h] am] edi 
pen p 


A condição de estacionaridade é 


sI= 22| | Ulm fm 9] (Loç)dacs tm 


E 


AS 


E | L( ComlfmTh] (Cond BS, | su 


BNB 


o desta condição resulta o sistema de equações lineares: 


(13) 


bà | | (Log) tl, jdáe+ | (C Fm) (Co,)d pa ii 
e lo ; 


Ae BNB 


uu 
lo 


(14) 


a s| Lo) gdaç+ | (Co )hdpo | 


A e BNB 
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Os dois métodos indicados podem considerar-se 
casos particulares do chamado método dos resíduos pe- 
sados [16], que utiliza um conjunto de funções w, e 
impõe que os integrais dos produtos do resíduo por cada 
uma delas se anulem, isto é, que 


» jut-gwda+/(Ct-hw dp=0 


Ú 


(15) 


“ 


A B 
Efectivamente, enquanto o método das colocações 
corresponde a utilizar funções w, que tomam o valor 


um num determinado nó de A e zero em todos os outros, 
o método dos minimos quadrados resulta de se fazer 


E Ph 
W 


IC, 


em à 


em 8 (16) 


como claramente mostra a equação (13). 

Um terceiro método, o método de Galerkine, con- 
sisto muito simplesmente em utilizar as próprias funções 
7 n como pesos, isto é, em fazer 


2| [iLt—-g)p,das+ [tet-meçdp'|=o 


[O 


A pen p 7 
Resulta o seguinte sistema de equações: 


E] - 


> | Lyn) pçdás + fc Pm! Pad | fm= 
hi a a | 
| cm 
= 0 o vqdact [od ar] 
e o d a á (18) 


Pat Ren B 


Outro tipo de métodos que apresentam vantagens 
grandes sobre os dos resíduos pesados é o dos que re- 
correm a princípios variacionais associados com a equa- 
ção diferencial (7) e respectiva condição de fronteira, 

Às vantagens compreendem-se muito bem se con- 
siderarmos um caso concreto que pode ser o do princi- 
pio da energia potencial total associado com as equações 
da Elasticidade. 

O princípio afirma que a solução exacta minimiza 
o funcional 


cm» 


F=U- [fu da-| udp 


A ; (18) 
PAN B, 
no conjunto dos campos compativeis. 
Ora, admitindo a lei de Hooke, sob a forma 
cj = Nx Ex (20) 
vem 
Far A 
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Como 
cu= 1/2 (u,;+ U;4) (22) 
[E] 
Vad Mix e AiS ETTO (23,24) 
temos ainda 
U= 122] hai, Ux, 16 À ii 
Á 


Uma das modalidades de utilização do método dos 
elementos finitos na Elasticidade (justamente a mais 
habitual) corresponde a empregar o princípio da energia 
potencial total e a discretizar os deslocamentos por in- 
termédio de 


4, Nato 126) 


Substituindo em (25) e, depois, em (19), obtém-se 


" 
o ” A é : 
Pr = 2 1 /2 | h tin] y ml hi mN.l da Um u KI 


c . 
ae 
ha] Em a = e 
=> Es | f. Fm O AC+ o; Ymd ç) Uim (27) 
r = a 
ar 5º nB, 
A condição de estacinaridade é 
cm 1 | = =» E ni 
SF = | 2) | Mia Pmj Pnd dd Ui 
ia 
Pai 
" à di | = : +, a 
ú 2 | | FP d+ | 0; 9md Bo ) 9 Um O 58) 
Ê a 4 a 
ae 5º 18, 
donde resulta, dada a arbitrariedade dos à Lim 
E temia Vin e Tim (29) 
em que 
K. =W | Naa Pai dAC 
tnimin cama nl ma “nn 
i Gado” (30) 
Ae 
fim > a! | ti vm d de + | a om 4º) 
e, id (31) 
a 8º ns, 
Observa-se que, no caso da isotropria [17]. 
q, = an Eu + AE É = [Za dk TA 0,9 «a E (32) 
e portanto, 
hj = 208,6) TAS; Sa (33) 


A equação (29), contrapartida discreta das equa- 
ções de Navier, não pode deixar de provocar uma obser- 
vação que julgamos oportuna e importante. 
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Essa observação é a de que, apesar de o formalis- 
mo matricial ter sido louvado durante tantos anos como 
o mais adequado à análise automática de estruturas, a 
experiência desses mesmos anos parece indicar que 
não é bem assim. 

A notação indicial tem com efeito vantagens apre- 
ciáveis sobre a matricial, não só norque são as opera- 
ções indiciais, e não as matriciais, que se programam 
directamente, mas também porque estas só podem subs- 
tituir as primeiras se não houver mais de dois índices 
afectos a cada grandeza. Ora, como se vê, até no caso 
presente, que está longe de competir em complexidade 
com aplicações mais sofisticadas do método dos elemen- 
tos finitos [18], aparecem grandezas com mais de dois 
indices. 

Por outro lado, a utilização de equações discretas, 
que talvez prenunciem uma Mecânica dos Meios Discre- 
tos formulável sem qualquer recurso a teorias continuas, 
substituem hoje o rude processo que consiste em progra- 
mar o cálculo da matriz de rigidez elemento a elemento, 
sem procurar uma lei de formação comum aos vários 
elementos, lei essa que o recurso ao formalismo indicial, 
tanto na Mecânica dos Meios Continuos como na técnica 
dos elementos finitos, torna normalmente fácil atingir. 

Mas voltemos às vantagens dos métodos variacio- 
nais. 

Verifica-se que a ordem máxima das derivadas 
das funções de forma que aparecem na expressão dos 
Kikmn é inferior à das que intervêm nas equações de 
Navier, isto é, nas equações diferenciais do problema 
que estamos a considerar. 

Isso significa que as funções de forma E - que 
são em geral polinominais e devem ter grau pelo menos 
igual ao das derivadas, podem ser mais simples. No 
caso presente, por exemplo, podem ser lineares, e não 
quadráticas, como deveriam ser se recorrêssemos aos 
métodos das colocações ou dos minimos quadrados. 

O método de Galerkine pode também adquirir esta 
vantagem se se transformar o primeiro membro de (18) 
com o auxílio do teorema da divergência. É de resto pos- 
sivel demonstrar que sempre que determinadas equações 
discretas forem obtidas por aplicação de um princípio 
variacional o poderão ser também por aplicação do 
método de Galerkine. 

O método de Galerkine é no entanto uma técnica 
mais geral que a dos principios variacionais, dado que 
estes nem sempra existem a o método de Galerkine é 
sempre aplicável. Sendo assim tão geral e gozando da 
vantagem quo indicâmos relativamente ao das coloca- 
ções, o dos minimos quadrados, compreende-se que estes 
últimos, aparentemente tão óbvios, nunca sejam apli- 
cados, 8 que o método dos elementos finitos seja quase 
sempre apresentado como técnica variacional quando, no 
fundo, ele não é mais que uma técnica de interpolação. 


2.4 — Comparação do método dos elementos finitos 
com o das diferenças finitas 
Estamos neste momento em condições de fazer 


entre os dois métodos uma comparação razoavelmente 
lúcida, 
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Examinemos primeiro a questão sob o ponto de 
vista da interpolação, 

Constatâmos a dificuldade em fazer interpolações 
globais sobre 4 com funções da classe C*. O método 
dos elementos finitos remove essa dificuldade sacrifi- 
cando a continuidade, o das diferenças finitas, sacrifi- 
cando a globalidade. Considera com efeito, para formar 
cada uma das equações do sistema, uma interpolação 
local, isto é, válida somente na vizinhança do nó a que 
a equação se refere, 

Observamos neste ponto que o método das dife- 
renças finitas não obriga a uma rede regular como muitas 
vezes se afirma. Pode na realidade fazer-se a interpola- 
ção em cada ponto por intermédio de uma função poli- 
nomial cujos coeficientes se calculam impondo a con- 
dição de tomar no próprio nó e nos nós vizinhos (a de- 
signar) os valores a interpolar. Esta mesma óptica mostra 
que as dificuldades habitualmente associadas à expres- 
são das condições de fronteira não são reais. 

Sob o ponto de vista da aproximação, o que pode- 
mos dizer é que os dois métodos não se distinguem. 
isto é, tanto os métodos dos resíduos pesados (coloca- 
ções, minimos quadrados, Galerkine, etc.), como os ba- 
seados na utilização de princípios variacionais, podem 
ser utilizados com a técnica das diferenças finitas, tal 
como com a dos elementos finitos. 

A única diferença reside em que tem sido costume 
utilizar a técnica das diferenças finitas com a das colo- 
cações e a dos elementos finitos com princípios varia- 
cionais ou com o método de Galerkine. 

Nada obsta no entanto a que se utilize a das dife- 
renças finitas também com o método de Galerkine ou 
com princípios variacionais. Cremos que este procedi- 
mento, que já tem sido recomendado por outros auto- 
res [19], constituirá passo essencial para a regeneração 
de uma técnica que tem sido imerecidamente caluniada. 

Observa-se somente, em abono dos elementos fi- 
nitos, que a localidade da interpolação associada com as 
diferenças finitas exige que o operador seja um operador 
diferencial, exigência essa que não aparece no caso dos 
elementos finitos. O forte sentido físico deste último 
método continuará por outro lado a favorecê-lo. 


2.5 — Exactidão contra generalidade 


Disse-se já que o preço da generalidade é a perda 
da exactidão. Tudo está pois em saber até que ponto o 
engenheiro pode dar-se ao luxo de nerder exactidão. 

Na prática, para responder a esta questão, o enge- 
nheiro tem que ter ideia do erro que comete ao abordar 
a resolução de um determinado tipo de problemas com o 
auxílio da técnica dos elementos finitos. 

Por outras palavras, o método não deve ser aplicado 
sem controle, e este ponto é tanto mais importante quanto 
ele se tem divulgado extraordinariamente e está hoje 
a ser aplicado por pessoas de reduzida experiência que 
consideram os resultados da análise estrutural de uma 
placa tão indiscutíveis como os da de um pórtico. 

Ora a verdade é que o método dos elementos fini- 
tos é daqueles instrumentos que não pode ser depositado 
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em todas as mãos, mas só naquelas que sirvam cabeças 
capazes de subdividir adequadamente o corpo em ele- 
mentos (e quanto a este ponto tem interesse referir as 
investigações tendentes a estabelecer a decomposição 
óptima [37]) e de adequadamente interpretar os resul- 
tados. 

Nisto, o método tem semelhanças com as técnicas 
experimentais, das quais tudo pode tirar-se se não hou- 
ver um mínimo de cuidado e de espírito crítico. E, como 
no caso destas, nenhum relatório apresentando um cál- 
culo efectuado com o auxílio do método dos elementos 
finitos terá valor se não indicar todas as condições 
em que foi executado, de modo, que possa ser repetido 
o criticado por qualquer leitor qualificado. É nomeada- 
mente indispensável que se indique como se fez a de- 
composição do corpo em elementos. 

Torna-se por outro lado muito importante comparar 
os resultados obtidos pelo método dos elementos finitos 
com resultados obtidos por métodos mais precisos, nos 
casos particulares em que estes possam ser aplicados, 
a extrapolar as conclusões para outros casos em que O 
primeiro seja o único viável. 

É bem típica deste procedimento a investigação 
[20] feita sob nossa orientação para controlar os resul- 
tados fornecidos pelo método dos elementos finitos para 
o problema da estaca cravada num terreno supostamente 
elástico e submetido à acção de uma força horizontal. 

A dificuldade do emprego do método dos elemen- 
tos finitos está neste caso associada com o facto de, 
sendo o domínio ilimitado, se correr o risco de utilizar 
um número insuficiente de elementos. 

O método por nós empregado, baseado na solução 
do problema de Mindlin, só pode servir no caso demasia- 
damente particular da homogeneidade e linearidade. Só 
para esse caso permite pois, em rigor, tirar conclusões 
sobre a conveniência de uma certa decomposição do 
semi-espaço em elementos. 

É porém natural admitir que uma mesma decompo- 
sição que sirva no caso linear e homogéneo possa ser 
utilizada na resolução de problemas com heterogeneida- 
des e não-linearidades, e até que o erro de discretização 
seja da mesma ordem de grandeza. 

Os métodos numéricos de alta precisão, como os 
métodos integrais, conservam pois a sua utilidade. 

Conservam-na também os métodos formais, capazes 
sobretudo de projectar luz sobre situações limites que os 
métodos numéricos só grosseiramente são susceptíveis 


de tratar. 
Quem, por exemplo, terá a coragem de concluir, de 


uma análise feita com elementos finitos, que num disco 
submetido à acção de forças concentradas iguais e opos- 
tas, actuando nas extremidades de um diâmetro, as ten- 
sões normais na direcção perpendicular ao das forças 
se distribuem uniformemente? E no entanto trata-se de 
conclusão importante, a avaliar pelo partido que dela 
tira o chamado «ensaio brasileiro [21]». 

As soluções formais devem pois continuar a ser 
utilizadas, isoladas ou sobrepostas com outras determi- 
nadas pelo método dos elementos finitos. Com sobrepo- 
sições deste tipo pode abordar-se noméadamente o es- 
tudo dos efeitos de fendas e de qualquer tipo de sin- 
gularidades [22]. 
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Deve pois concluir-se que, apesar de todas as suas 
qualidades, o método dos elementos finitos não é uma 
panaceia, e que não se justifica de modo nenhum a 
atitude daqueles que, desde que ele existe, julgam possi- 
vel desprezar a análise. 


3 — PARA ALÉM DA DISCRETIZAÇÃO 
3.1 — Investigações em curso 


Quem se der ao trabalho da examinar o índice de 
uma revista como o «International Journal of Computer 
and Structures» terá a surpresa de constatar que muito 
poucos artigos têm por tema o método dos elementos 
finitos, o que de resto explica que, embora as investiga- 
ções consagradas aos métodos dos elementos finitos 
tenham vindo a diminuir de volume, o número de traba- 
lhos publicados sobre a análise e síntese de estruturas 
mostra forte tendência para aumentar [23]. 

O método dos elementos finitos não é na realidade 
senão uma técnica de discretização e, portanto, nada 
mais que um instrumento matemático para o estudo de 
fenómenos complexos de cuja simulação e consequente 
análise se ocupam hoja em dia cada vez mais os inves- 
tigadores. 

Temos assim, além das investigações dedicadas 
ao método dos elementos finitos, trabalhos que se dis- 
tribuem nomeadamente pelos seguintes domínios: 


—— Qutros métodos numéricos; 
— Idealização e análise de estruturas especiais; 


—— Teorias elásticas continuas lineares (vigas, cas- 
cas, maciços, etc.); 


— Termoelasticidade; 

—. Elasticidade finita; 

— Problemas geometricamente não-lineares; 
— Não-linearidades físicas; 
—. Plasticidade; 

— Viscosidade; 

—. Instabilidade; 

— Dinâmica; 

—. Solicitações; 

—. Segurança; 

— Optimização; 

— Programação. 


Não se incluem aqui os estudos feitos sobre mate- 
riais, em relação aos quais há sobretudo que mencionar, 
por um lado, as contribuições dos engenheiros, e, no 
caso que nos interessa especialmente, que é o da enge- 
nharia civil, os estudos efectuados sobre o material be- 
tão armado, por outro lado, as contribuições dos físicos 
8 dos matemáticos que deram considerável desenvolvi- 
mento à moderna Mecânica dos Meios Contínuos. 

Na evidente impossibilidade de tocar em todos os 
domínios mencionados, limitar-nos-emos a falar um 
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pouco sobre os problemas geometricamente não-lineares, 
porque julgamos poder comunicar algo não inteiramente 
banal, e sobra a optimização, dando informação que, 
para o capítulo seguinte, será de utilidade. 


3.2 — Problemas geometricamente não lineares 


A equação (29), baseada na equação (22), não é 
válida senão enquanto forem admissíveis as hipóteses 
simplificativas que caem dentro daquilo a que se cos- 
tuma chamar linearidade geométrica. 

Em casos em que não sejam de admitir tais hipó- 
teses, haverá que partir da expressão exacta das rela- 
ções deformações — deslocamentos, isto é, 


Ej = 1/2 (u;;+ U;; + U cj Ui ) (34) 


Ria 


que, introduzida na expressão da energia potencial total 
(19) juntamente com (26), conduz à seguinte condi- 
ção de estacionaridade [24]: 


| Aimnttom + Pen Upp) P snPnl JR ii 


A E 


NE 


| [+ f ' : HF | -— 4 
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a 
(35) 


A equação (35) substitui a equação (29) tanto nos 
casos em que intervenham grandes deslocamentos como 
grandes deformações, 

Em trabalhos anteriores [25] preconizámos no en- 
tanto um método apropriado para grandes deslocamentos 
8 pequenas deformações, que é na realidade muito mais 
simples que o baseado na resolução da equação (35). 

Esse método considera os deslocamentos Em de- 
compostos em duas parcelas, representando a primeira 
um deslocamento de corpo rígido de componentes à in 
tal que sejam pequenos os deslocamentos da segunda 


parcela, u - Tem-se pois 


KIM 


Up = Ur + Un (36) 


n 


Sendo os deslocamentos Un 
utilizar a equação (29) sob a forma 


pequenos, podemos 


mm 


Kikmn Un = f, (37) 


uma 


Introduzindo (36) em (47), obtém-se 


se f 
Kim Un "im Ed im (38) 


em que 


É «E ” 
f imo Kikmn Un (39) 


Quer dizer: a equação (29) pode ser usada desde 
que haja o cuidado de somar às forças reais uma parcela 


de forças fictícias e! 


im dependente do deslocamento de 
compo rígido 0... 
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A dificuldade em utilizar a equação (38) reside 
evidentemente na determinação dos deslocamentos à a 

A dificuldade não chega a existir se o problema 
consistir em calcular as forças a partir de dada configu- 
ração. Se o problema a resolver for no entanto o inverso, 
o de calcular os deslocamentos a partir das forças, po- 
der-se-á recorrer a um processo incremental em que os 
deslocamentos rígidos e os coeficientes de rigidez K ,.un 
correspondentes a cada configuração poderão ser deter- 
minados a partir da configuração anterior. 


3.3 — Optimização 


Parece dever-se a Galileu o primeiro trabalho sobre 
optimização estrutural, datado de 1638 e incidindo so- 
bre vigas de igual resistência à flexão [26]. Seguiu-se- 
-lha uma quinzena de autores, até 1904, ano em que 
Mitchell publicou o seu notável trabalho, «The Limit of 
Economy of Material in Frame Structures» [27], consa- 
grado à determinação da geometria óntima de treliças. 

Talvez nenhum engenheiro tivesse tido a ideia de 
projectar estruturas semelhantes às de Mitchell, até por- 
quo a teoria negligencia os fenómenos de instabilidade 
que condicionam fortemente as barras comprimidas. Não 
deixaremos no entanto de observar que certas estruturas 
biológicas se parecem notavelmente com estruturas de 
Mitchell. 

Os processos actuais de optimização de estruturas 
podem ser classificados nas seguintes categorias [28]: 

a) optimização geométrica; 

b) optimização baseada em critérios de optimali- 

lidade; 

c) optimização por via numérica; 

d) métodos mistos. 


Cabe na primeira a já mencionada teoria de Mi- 
tchell. 

Os critérios de optimalidade são, em princípio, 
condições a que a solução óntima tem que obedecer, 
isto é, condições necessárias de óptimo. 

Às vezes, no entanto, tais critérios não são mais 
que condições que caracterizam, não a solução óptima, 
mas soluções que se sabe estarem suficientemente pró- 
ximas da óptima. 

No que se refere à optimização estrutural, o exemplo 
clássico, no caso de a função objectivo ser o peso da 
estrutura ou de se fixar a geometria desta, é o critério 
da simultaniedade dos modos de colapso, o qual supõe 
que a estrutura óptima é aquela para a qual todos os 
possíveis modos de colapso se atingem simultanea- 
mente [29]. 

Se o número de modos de colapso for igual ao de 
variáveis definidoras da estrutura, o critério fornece um 
sistema com o mesmo número de equações e de incóg- 
nitas cuja solução determina a estrutura óntima. Se for 
inferior, não chega evidentemente para determinar a es- 
trutura óptima, mas facilita a sua determinação vor via 
numérica, servindo assim de base a um método misto, 


O critério do aproveitamento máximo das secções 
(«fully-stressed design»), implícito na investigação de 
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formas de igual resistência, e que no caso geral de sis- 
temas da elementos finitos se transforma em critério do 
aproveitamento máximo dos elementos, constitui parti- 
cularização do anterior [30]. 

Trata-se de um critério cuja legitimidade é evidente 
para as estruturas isostáticas, em que os esforços nas 
barras não dependem das respectivas secções, mas que 
deixa de ser válido para estruturas hinerstáticas em que 
a solução de peso mínimo nem sempre é uma solução 
com aproveitamento máximo das secções. Um teorema 
demonstrado em 1873 por Maurice Lévy afirma mesmo 
a inexistência em geral, isto é, para solicitações arbitrá- 
rias, de treliças hiperstáticas de igual resistência [261]. 

A experiência mostra no entanto que, pondo de lado 
não só as treliças como certos casos patológicos, e desde 
que não se imponham outros contrangimentos além do das 
tensões máximas, o «fully-stressed design» dá estrutu- 
ras suficientemente próximas das óptimas. Grandes van- 
tagens do critério, que é normalmente utilizado como 
base do um processo iterativo, são a rapidez de con- 
vergência do referido processo, rapidez essa que se veri- 
fica sempre que se estiver fora dos tais casos patoló- 
gicos, e o facto de o número de iterações necessárias 
não depender praticamente do número de variáveis em 
jogo. 

A ideia do aproveitamento máximo é no fundo a 
de que, na solução óptima, cada uma das desigualdades 
que traduzem os constrangimentos das tensões devem 
ser satisfeitas como igualdades. Não admira que tal 
concepção tenha sido generalizada e aplicada a outros 
tipos de constrangimentos, como os que se traduzem por 
limitações dos deslocamentos, 

Tem-se apresentado também critérios de optimali- 
dade baseados nos teoremas de mínimo da Teoria das 
Estruturas. Principalmente devidos a Prager 31, 32. 
33], consideram como única função objectivo o peso 
(ou o volume) da estrutura, mas contemplam diferentes 
tipos de constrangimentos. 

Todos os critários mencionados supõem fixa a geo- 
metria das estruturas. 

Os problemas que envolvem a determinação de 
geometrias óptimas («optimal layout») [34] são muito 
mais complicados, não se conhecendo mesmo para eles 
critérios gerais de optimalidade. 

A cptimização numérica [35], quase sempre ba- 
seada em técnicas de programação matemática, linear, 
quadrática, dinâmica ou geométrica, é então em regra 
a única possível, dado que a optimização geométrica, 
baseada na teoria de Mitchell, é de alcance muito 
limitado. 

Mesmo assim o problema da determinação da con- 
figuração óptima só foi resolvido até agora, conside- 
rando ou não o risco da instabilidade, para treliças estati- 
camente determinadas. Demonstrou-se por outro lado 
que, sendo as tensões máximas independentes dos es- 
forços, o que só acontece quando não se considera a 
instabilidade, a solução existe sempre mas não é ne- 
cessariamente única. 

Os resultados que acabam de ser mencionados 
supõem uma só hipótese de carga e constrangimentos 
unicamente nas tensões. Em casos em que não se déem 
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tais condições, prova-se que a solução óptima é uma 
traliça hiperstática de barras não maximamente aprovei- 
tadas, o que de resto concorda com o já citado teorema 
de Maurice Lévy, 

Menciona-se como curiosidade uma técnica recente 
para a determinação da configuração óptima de treliças, 
a qual se baseia no teorema de Mitchell, que fornece 
a configuração inicial, e nos princípios biológicos da evo- 
lução, que sugerem estratégias óptimas para processos 
da optimização de sistemas [36]. 

Aparenta-se com este problema da determinação 
de configurações óptimas o problema, já mencionado no 
capitulo anterior, da determinação da decomposição 
óptima da um corpo em elementos finitos, para o qual 
foi já estabelecido um critério de optimalidade que con- 
sisto em dispor os elementos de harmonia com as linhas 
(ou, no espaço, superfícies) de igual densidade de ener- 
gia de deformação [37]. O problema pode também 
evidentemente ser abordado exclusivamente por técnicas 
numéricas. 

Os métodos numéricos têm em todos estes proble- 
mas dae optimização a grande vantagem da generalidade, 
generalidade essa que lhes permite considerar funções- 
-objectivo não mecânicas, em relação com as quais não 
poderão já estabelecer-se critérios de optimalidade 
baseados nos princípios da Mecânica. A necessidade 
de considerar tais funções-objectivo aparece sobretudo 
na engenharia civil, visto que na aeronáutica a minimi- 
zação do peso da estrutura tem interesse bem real, 

De qualquer modo, o que parece lógico é encarar 
o problema da segurança como problema de optimiza- 
ção [38], tomando para função-objectivo o custo gene- 
ralizado ou total, 


c=c+pe. (40) 


resultanta da soma do custo inicial, cj... do produto 
da probabilidade de colapso, p, pelo custo do colapso, 
C. + por sua vez resultante do custo dos prejuizos, que 
á independente de p, somado com o custo inicial, isto é, 


= C+ec, (41) 


ia 


Observa-se que o custo inicial depende de pn e da 
configuração E da estrutura, pelo que a função-objectivo 
podo ser expressa por 


clppE)-=c Ip E) (1+p)+p e, (42) 


O problema está pois todo na determinação da 
configuração e da probabilidade de colapso que dêem 
origem ao custo generalizado mínimo. Observa-se que 
a configuração é também função de op, dado que a pro- 
babilidade de colapso condiciona os coeficientes de 
segurança, 

Em regra, para simplificar o problema, fixa-se a 
probabilidade de colapso, p, e determina-se a configu- 
ração correspondente ao mínimo custo generalizado. 
Não chega portanto a determinar-se o valor de p que 
conduz ao valor minimo dos mínimos de c. Outras vezes 
ainda, fixa-se o custo generalizado e determina-se a 
configuração que conduz ao valor mínimo de »p. 

De qualquer modo, quanto mais realisticamente se 
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põe o problema da optimização das estruturas, mais 
formidáveis são as dificuldades que se levantam no 
caminho da sua resolução. Ao tentar libertar-se da fama 
mais ou menos merecida de somente resolver problemas 
académicos, a síntese estrutural vê assim surgir obstá- 
culos dos que normalmente não se removem sem uma 
completa reformulação dos problemas, neste caso só 
possivel depois de uma mais profunda reflexão sobre 
a utilização dos computadores digitais na Engenharia. 
A ela so dedica o próximo capítulo, 


4 — PARA ALEM DOS CALCULOS 
41 — O papel do engenheiro e do computador 


A característica essencial do engenheiro é a de to- 
mar decisões, relativas por exemplo à localização das 
obras, ao seu dimensionamento, aos materiais, às des- 
pesas, à segurança [39]. 

Para tomar tais decisões, o engenheiro baseia-se 
em informações relativas a exigências mínimas a que 
as soluções devem satisfazer, a códigos a que é neces- 
sário obedecer, a condicionamentos de ordem topográ- 
fica o geológica, aos recursos disponíveis e a limitações 
próprias dos materiais, 

Um gabinete de projectos pode justamente ser 
visto como um sistema de processamento da informa- 
ção, onde a informação bruta se condensa, transforma, 
cataloga, memoriza e combina com mais informação. Tal 
processamento ou manipulação da informação faz-se 
através do seu arquivo, classificação e comunicação, 
da execução de cálculos, da verificação e ensaio de so- 
luções, da comparação com soluções antigas, do estudo 
de variantes, da elaboração de desenhos, gráficos e 
tabelas, 

É nestas tarefas que a maior parte dos engenheiros 
exerco actualmente a sua actividade. Quer dizer: a maior 
parto dos engenheiros são hoje em dia meros processa- 
dores da informação desempenhando funções que seriam 
com absoluta vantagem desempenhadas Dor computado- 
res, 

No futuro, o engenheiro aparecerá cada vez mais 
como agente de decisões e o computador como proces- 
sador da informação em que vão basear-se essas deci- 
sões, 


4.2 — Os problemas de engenharia como problema de 
optimização 


Os problemas do engenheiro são no fundo proble- 
mas de optimização [39]. 

Ao engenheiro não interessa com efeito somente 
uma solução, mas a solução óptima. À ideia de optimiza- 
ção só não chega a ser explicitada nos casos mais cor- 
rentes porque uma longa experiência leva a adoptar direc- 
tamente soluções que se crê implicitamente estarem 
próximas das óptimas. 

Duas metodologias são possíveis. 

A primeira consiste em formular o problema em 
termos matemáticos precisos, resolvendo-o por técnicas 
matemáticas conhecidas, como as qua se mencionaram 
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no capítulo anterior, tudo se passando dentro do com- 
putador sem intervenção humana desde que os dados 
entram até que os resultados saem (optimização di- 
recta). 

Este método tem dois inconvenientes graves: a di- 
ficuldade da adaptação dos problemas reais a processos 
sistemáticos de optimização, e a dificuldade da própria 
formulação matemática dos problemas, mais precisa- 
mente, da fixação das funções objectivo, isto é, daquilo 
que se pretende optimizar. 

Muitas vezes, com efeito, a solução resulta de um 
compromisso entre, por exemplo, a optimização econó- 
mica e a estética, ou entre o custo e a segurança. Ora, 
quase sempre, não só é impossivel formalizar esse com- 
promisso como cada uma das optimizações parciais é 
discutível em si mesma. Que significa optimização es- 
tética? Como poderão avaliar-se exactamente os riscos? 
Como poderão encerrar-se numa fórmula todas as pos- 
síveis implicações de determinadas decisões? 

A segunda metodologia (optimização por simulação) 
tem também na base, como é natural, um modelo ma- 
temático que simula o funcionamento do sistema a opti- 
mizar. Mas, enquanto na optimização directa o modelo 
matemático fica subordinado à formulação e resolução 
do problema da optimização, nesta, basta apenas simular 
com a maior fidelidade e com suficiente generalidade 
par que todos os parâmetros significativos possam ser 
variados, As decisões sobre a variação dos parâmetros 
e o próprio reconhecimento da decisão óptima resultam 
da intervenção directa do engenheiro. 

Asim, neste último método, a sensibilidade e a ex- 
periência pessoal, por um lado, a tradição profissional 
por outro, permitem curto-circuitar o processo de opti- 
mização, levando a soluções que se admite estarem pró- 
ximas das óptimas. 

Estarão? 


Se nalguns casos isso é duvidoso, noutros a sensi- 
bilidade diz-nos que efectivamente estão próximas dum 
óptimo que por resultar de um compromisso, às vezes 
nebuloso, nem sempre somos capazes de definir com 
precisão. 

A conclusão importante a reter é a de que a auto- 
matização completa não é geralmente possível em proble- 
mas reais, embora o possa ser em problemas parcelares, 
a utilização da cujas soluções fica no entanto sempre 
dependente de um julgamento lúcido por parte do enge- 
nheiro. 

A optimização por simulação exige no entanto uma 
fácil comunicação homem-máquina [40], o que envolve 
o desenvolvimento da linguagens de programação ade- 
quadas, fácil acesso aos computadores («time-sharing», 
etc.) e dispositivos adequados para entrada de dados 
& saída de resultados, como pantalhas de raios catódi- 
cos, traçadores de gráficos, dispositivos de leitura de 
desenhos, máquinas de desenhar, etc. E há certamente 
mais a esperar, neste momento, do desenvolvimento 
acertado e da divulgação de dispositivos deste tipo que 
do progresso dos processadores centrais que costumam 
ser considerados como a parte nobre dos computadores. 


4.3 —- A adaptação da profissão aos computadores 
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Embora o método clássico do engenheiro apresente, 
como se viu, caracteristicas que não podem ser sacrifi- 
cadas, não há dúvida de que a profissão precisa de ser 
adaptada às circunstâncias novas e, muito especial- 
mente, aos computadores. 

Uma primeira adaptação diz respeito aos gabinetes 
de projectos e às empresas em geral. 

Um gabinete de projectos tem de passar a consti- 
tuir um sistema integrado, como um navio, e o tipo de 
organização antiga em que coexistiam, mas não comuni- 
cavam, um sector técnico e um sector comercial, depen- 
dentes ambos de uma direcção a cujo nível se combi- 
navam as informações que ambos forneciam, tende a 
ser substituído por um modelo com um sector técnico- 
-comercial, orientado para os projectos, um sector ana- 
lítico, orientado para os problemas, e um sector de cál- 
culo, contendo o computador e os programadores, todos 
trocando informações entre si e com a direcção, que 
não deixa evidentemente de estar no centro, animando, 
coordenando e sintetizando [41]. 

Esta organização moderna tem evidentemente in- 
fluência sobre a formação dos engenheiros, suscitando 
O aparecimento de um novo tipo de profissionais, os en- 
genheiros analistas, adequados ao sector analítico e por- 
tanto voltados para as formulações gerais. Os especia- 
listas de tipo clássico,que estabelecem a ponte para os 
problemas concretos, actuam no sector técnico-comer- 
cial. Do contacto entre uns e outros nascem os modelos 
físicos dos problemas concretos, assim como do con- 
tacto entre o sector analítico e o centro de cálculo nas- 
cem os modelos matemáticos imediatamente programá- 
veis. É evidente que nas organizações mais rudimentares 
as funções serão acumuladas e não chegarão forçosa- 
mente a explicitar-se em sectores separados, 

Uma segunda adaptação diz respeito à codificação 
e formalização dos processos e à coordenação dos 
dados. 

Exige-se com efeito uma minuciosa reflexão sobre 
os processos tradicionais que permita simplificar, codi- 
ficar e formalizar, de modo a melhorar a comunicação 
entre empresas e a facilitar o acesso a dados sobre ma- 
teriais, produtos, regulamentos, normas, custos, etc [42]. 

Há portanto que estabelecer esquemas normaliza- 
dos para a execução de tarefas dent:o do processo de 
planeamento da const ução, passando pelo projecto, pela 
adjudicação e pela construção propriamente dita, poden- 
do tais esquemas ser eventualmente ap.esentados sob a 
forma da organogramas [43]. 

A utilização e aperfeiçoamento de sistemas de codi- 
ficação permitirá por outro lado que os desenhos conte- 
nham toda a informação [44] necessária e que ao 
mesmo tempo se mantenham simples, libertando o pro- 
jectista da tarefa de pormenorizar. O ideal será que ve- 
nham a ser lidos directamente pelos computadores. 


44 — À valorização de tarefas tradicionalmente 
menosprezadas 


Em virtude do que se tem dito torna-se claro que a 
revolução que os computadores vieram provocar tem van- 
tagens com significado em si mesmas, isto é, que seriam 
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válidas mesmo seg;mnão houvesse computadores A mi- 
núciosa reflexãboqueio aparecimento dos computadores 
determinou veio com efeito revelarcuma necessidade de 
racionalização que há: muito “tempo se fazia sentir e que 
esso aparecimento simplesmente tórnou mais: premente: 
orássimo se há calguns anos se tivesse: pergun- 
tado aos engenheiros qual cera “o produto mais im- 
portante: dos: gabinetes: de projecto; a maioria: talvez 
respondesse que eram os  cálculos,-& esta resposta fa- 
cilmente se explica se secatender a rquerera nos cálculos 
que essds: np sd a maior parta: do seu 
tempos su lovin o 5 di eu Gl dr na 
Agora, porém, quera traa: são anbdeár cálculos; 
como processamentorde informação, tende'a passar aos 
computadores, começa-se: a reconhecer que os produtos 
acabados dos: gabinetes de! projectos são essencialmente 
os desenhos, os: cadernos-de encargos e os orçamentos. 
Esta- conclusão tem qualquer coisa de irónico por- 
que; se os desenhos sempre foram considerados impor= 
tantes, a execução de |mbdições e'o cálculo de custos 
unitários para a elaboração de orçamentos ve a 'prepara- 
ção de cadernos: de encargós eramrconsiderados tarefas 
pouco nobres. | Do | pu sb amiga 
Encarar estas tarefas na sua generalidade, de modo 
a equacionar o problema da sua automatização, teve no 
entanto: a consequência de as nobilitar. | 
roriTudo se cprepara-hoje para que o computador possa: 
ler automaticamente os desenhos codificados, formando 
artista: dos trabalhos a efectuar e dos: números que tra- 
duzem as: respectivas: medições: Recorrendo aos custos 
unitários: conservados em memória: [45], 'o computador 
elabora então, por um: lado; automaticamente vos“ orça- 
mentos: Por outro lado, a meésmartista de-tarefas permi+ 
te-lhe reunir -e! imprimir um ocerto número de disposi- 
sa E id [46] e formar 'o- caderno de iii 
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4.5-— À valorização da intuição ce da imaginação 

“Sejaragora O nosso último bonto de reflexão. | 

. Automatizando, o computador liberta O engenheiro 
das tarefas dé rotina, pêrmitindo-lhe' Pensar em termos 
globais)" ter uma visão de conjunto, educar: E] sensibili- 
dado mediânte a exploração fácil, para cada modelo que 
simúle à réatidade, do espaço das soluções possíveis. 

“Nestas condições, à construtor | pode libertar-se ver 
dadeiramente, como atrás 'se afirmou, “do segundo incon- 
veniente trazido pela primeira revolução: : a - decadência 
da intúlição” 5 dá Di da à € OB9ESUUM À 
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Máquina de turing para converter expressões 


em notação “Infix” em expressões 


em notação “Postfix” (<) 


1. Definição da máquina e redução de estados 


RESUMO 


Este artigo é o primeiro de dois dedicados ao pro- 
jecto de uma máquina de Turing para converter notação 
intix em postfix. São amplamente utilizados, como auxi- 
liares no projecto, métodos de análise de estrutura al- 
gébrica de máquinas sequenciais, que têm vindo a ser 
estabelecidos no dominio da Teoria dos Automata. 

Neste primeiro artigo é apresentado um algoritmo, 
definida uma máquina de Turing nele baseada, e redu- 
zidos os estados da máquina de modo a obter uma má- 
quina minima. 


Tintrodução 


Existe um grande conjunto de resultados relativos 
à estrutura algébrica de máquinas sequenciais, recente- 
mente estabelecidos no dominio da Teoria dos Automata. 
Estes resultados nem sempre são utilizados em projectos 
de máquinas sequenciais. À principal finalidade deste 
trabalho é ilustrar, através de um exemplo, a utilidade 
que pode ter, num projecto deste tipo, a aplicação de 
métodos de análise da estrutura algébrica de máquinas. 

O exemplo em discussão será uma máquina de Tu- 
ring para converter expressões algébricas de notação 
infix (a notação usual) para notação postfix (ou polaca). 
Na notação postfix as operações são indicadas por or- 
dem de execução e depois dos operandos respectivos; 
assim, a conversão infix-postfix é realizada em compila- 
dores para processamento de expressões algébricas e, de 
forma implícita, na utilização de certas máquinas de 
calcular. 

Parece-nos conveniente começar por sistematizar 
O que se entenderá por notação infix e por notação post- 
fix. 


L. T. MAGALHAES 
M. P. O, RICOU 


ABSTRACT 


This is the first of two papers related to the design 
oí a Turing machine for the translation of infix notation. 
to postfix notation. Methods for the analysis of algebraic 
structure of sequential machines, recently developed as 
part of Automata Theory, are largely used. 

In this first paper we present an algorithm, define 
a Turing machine realizing it, and reduce machine states 
in order to get a minimal machine. 


2 Notação infix 


1. As expressões consistem numa sucessão de ope- 
randos, que podem ser precedidos por operações uná- 
rias, tendo entre si um símbolo de operação binária. São 
iniciadas e terminadas pelo símbolo =. 


2. Os operandos são caracteres alfanuméricos 
isolados, exceptuando-se o caracter F, ou o resultado de 
operações (binárias ou unárias) realizadas sobre ope- 
randos. 

3. As operações binárias são: 

a) adição + 

b) subtracção ni 

c) multiplicação ” 

d) divisão f : exemplo: A/B «+ A : B 
e) exponenciação + : exemplo: AtB-> AB 


4. As operações unárias são representadas por 
uma sequência de 4 caracteres alfanuméricos, dos quais 


(*) Este trabalho foi realizado em 1973/74 enquanto alunos da cadeira Sistemas Digitais — IST, Lisboa, 
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o primeiro é sempre F. O operando a que dizem respeito 
é o que está imediatamente a seguir e deve estar entre 
parêntesis curvo, 


Exemplos: FSIN(X), FEXP(A“(B/C)) 


5. As operações binárias actuam sobre os dois 
operandos entre os quais se encontram, e na ordem em 
que figuram na expressão. Consideram-se certas regras 
para a ordem de operações: 


a) a ordem de prioridade das operações é 
ij + li) F,/ il) +, — 


b) quando uma subsucessão de operandos é ligada 
por operações com a mesma prioridade, a ordem 
de operação é do princípio para o fim. 


c) a ordem de operação assim estabelecida pode 
ser alterada pelos simbolos de parêntesis curvo 
(ou); as operações entre parêntesis devem ser 
executadas antes das operações que precedem 
e/ou seguem os parêntesis correspondentes. 


3. Notação postfix 


1. As expressões consistem numa sucessão de sim- 
bolos alfanuméricos isolados, operações unárias e biná- 
rias, São iniciadas e terminadas pelo simbolo =. 

2 
como para notação infix 


4. As operações unárias são representadas por uma 
sequência de 4 caracteres alfanuméricos, dos quais o pri- 
meiro deve ser F. O operando a que dizem respeito é o 
que está imediatamente antes. 


5. As operações binárias são consideradas como 
operando sobre os dois operandos que imediatamente 
as precedem e na ordem em que figuram na expressão. 


6. A ordem de operação é sempre do princípio para 
o fim da expressão, 


Nota: não exista uma ordem de prioridade de ope- 
rações, nem utilização de parêntesis. 


4. Exemplo de tradução de expressões 


A expressão algébrica 
= 2*3 1 (A+ B) + FSIN(FLOGIX)) + 2 = 
é traduzida para 


= 23AB + + *XFLOGFSIN + 2 + = 


5. Definição de máquina de Turing 
A definição formal adoptada para máquina de Tu- 


ring não é a mesma nara todosos autores [1],[2],[4]. 
Adoptaremos a seguinte definição: 
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— Uma máquina de Turing com n fitas é um sex- 
tuplo ordenado 


T=[(1,S,0,5,A, qo) tal que: 


1) 1= A” onde À é um conjunto não vazio e fi- 
nito de símbolos. | é o «conjunto das entradas». 


2) S é um conjunto não vazio e finito; o «conjunto 
dos estados». 


3) O=(WUAUÍI])xX le d,e,hj)n é o «con 
junto das saídas». d e e correspondem aos co- 
mandos para posicionar a cabeça de leitura na 
posição à direita ou à esquerda (respectiva- 
mente) da posição ocupada, é corresponde ao 
comando para manter a cabeça imóvel, e h cor- 
responde ao comando de paragem quando a 
operação está completa. | é o «símbolo de 
indiferença». Ei 


4) 5 é uma aplicação de um subconjunto não va- 
zio de |X S em S; a «aplicação estado-se- 
guinte». 


5) à é uma aplicação de um subconjunto não va- 
zio doe | X S em O: a «aplicação saida». 


6) q, e S é o estado inicial da máquina. 


6. Discussão de algumas alternativas 


A ordem das operações na expressão postfix à 
saida, depende dos critérios de prioridade entre as ope- 
rações e das alterações resultantes da inclusão de pa- 
rêntesis na expressão infix à entrada; a ordem dos ope- 
randos alfanuméricos é a mesma em ambas expressões. 

A memória necessária para decidir a ordem das 
operações na saída é reduzida se fôr utilizada uma zona 
de registo auxiliar, onde sejam memorizadas as opera- 
ções já lidas na expressão de entrada (pela ordem em 
que foram lidas a com a inclusão de parêntesis) e ainda 
não escritas na expressão de saida. A inclusão desta 
zona auxiliar permite reduzir o número de movimentos 
da cabeça sobre a expressão de entrada (e portanto o 
número de estados da máquina) e o número de simbolos 
que é necessário introduzir para aumentar o alfabeto 
da máquina. 

Optando pela inclusão desta zona de registo auxi- 
liar, é necessário considerar três zonas de registo: 


1) para a expressão de entrada 
2) para o registo auxiliar 


3) para a expressão de saída 


É simples verificar que a máquina pode ser reali- 
zada apenas com uma fita. Pode-se mesmo acrescentar 
que, neste caso, é possível projectar uma máquina de 
modo a que a zona da fita ocunada em qualquer ins- 
tante do processamento seja menor ou igual à zona ocu- 
pada pela expressão inicial. De facto, em qualquer ins- 
tante do processamento, as zonas 2) e 3) contêm um 
número de símbolos menor ou igual ao número de sim- 
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bolos já lidos na expressão de entrada; assim, as zonas 
2) e 3) podem sobrepôr-se à parte da zona 1) já lida. 

À utilização de uma única fita conduz a uma fre- 
quento posicionação da cabeça em cada uma das três 
zonas de registo e a uma memorização dos simbolos 
lidos durante os deslocamentos, com o consequente au- 
mento do número de estados. 

Visto que convém reduzir o número de estados 
[sem aumentar o número de símbolos do alfabeto), de 
modo a facilitar a atribuição de variáveis de estado e a 
codificação dos símbolos, projectar-seá uma máquina 
com três fitas (uma para cada zona de registo). 


7. Algoritmo para conversão infix-postlix 


Apresentamos na fig. 1 um algoritmo na forma de 
fluxograma. 


Convenções: 


1) x-»y significa que o símbolo lido na fita x é 
escrito na fita y. | 


2) xd, xe significam que a cabeça da fita x se 
move para a direita ou esquerda, respectiva- 
mente. 


3) a é utilizado para designar qualquer alfanumé- 
Fico. 


4) quando se não indicam transferências de sím- 
bolos como. convencionado em 1), entende-se 
quo não há alteração de simbolos; quando se 
não indica o movimento de uma cabeça de fita 
como convencionado em 2), entende-se que 
essa cabeça se não move. 


B. Diagrama de estados da máquina 


Com base no algoritmo apresentado é possível 
construir um diagrama de estados para a máquina (fig.2). 


Convenções: 

1) X | Aa representa a entrada (X, Y, Z) e a 
Y | Bb saida (A, a, B, b, €C, cj). 
Z | CE 

2) O representa uma qualquer operação. 

3) X representa o complementar de X. 


9. Máquina de Turing correspondente ao diagrama de 
estados apresentado 


A máquina representada pelo diagrama de estados 
anterior é uma máquina de Turing com 3 fitas T= (|,5, 
O, 5, à, q,), onde: 


— O conjunto 4, referido na definição (8 5), é o 
conjunto constituido pelos caracteres alfanumé- 
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ricos (36), pelos símbolos das operações biná- 
rias (5), e pelos símbolos (,), =, prefazendo 
um total de 44 simbolos. 


— S =. [Qu Qi. E: Gu) 


Falta ainda definir as aplicações 5 e à. Note-se que 
| tem (44)' elementos. Assim, se à e À fossem definidas 
extensivamento por tabelas, teriam de ter dimensão 8 X 
x (44)'. Felizmente é possivel ultrapassar esta difi- 
culdade. 


Sendo Op = [+,—,*,/,1), e « o conjunto dos 
símbolos alfanuméricos, considere-se a partição de |: 


« Es), definida por: 


; aa 
Ê — | Eu gds. 


= (=.(.27): ze A) 

([l=,=,2): ze À) 

ll=,y.2): yeA- (=) 0, ze A) 

[Mt,t,z): ze À) 

[ly z): yeOp— (1), ze A; 

(ley z): xe(+t,—), ye Op, ze A) 
=[(x,y,2): xel*,/), vyeOp— |+,—), ze A) 
= Iflxya xe(fº./h ye ld;—), ZeÃA) 

o= ([(x,y,27): x2eOp, yea, Z,2 A) 


PE Frio art da dd 
| E ha [P ba 
pum mM 


= 


En= [l(x,=,27): xe Op, z2 À) 


En= (lx, (,2):xe Op, ze A) 
En=(()l,2): ze À) 
Eu=MD)yz):vyea z e A) 


veOp Ul=), ze A) 
yeA—()), ze A) 

Xew yea, Ze À) 

xea, vyeOp Ul=,(), ze A) 
xeÃ, z cA) 


Em = ()y,z): 
Em= Ml, y,z): 
Es = (x,y, 2): 
Ev = [lx y,2z): 
Es = |(x,), 2): 


Note-se que, se q(0) e «w(0) designarem as par- 
tições nulas de S e de O, (=, q (0)) é um par de par- 
tição |-S [3] e (*, o(0)) é um par de partição |-0, como 
se pode concluir analizando o diagrama de estados apre- 
sentado. 

Podemos então pensar em definir duas funções: 


5*, de um subconjunto de E XS, e A*, de um sub- 
conjunto de: X S em O, de tal modo que 


8* (5 .q))=8 (q; |) sse xe E 


Ny a)=Atxa) sse Xe E 


Analizando o diagrama de estados apresentado foi 
possível chegar às tabelas | e Il, respectivamente para 
5* e A*, Atendendo às relações entre 5* e 5* e entre 
X* e X, ficam assim definidas à e À. 

O domínio de 8* e de Nº é 0 conjunto de elemen- 
tos de E X S para os quais está definido um valor nas 
tabelas. As entradas não preenchidas correspondem a 
situações de indiferença (don't care). 

Esta máquina foi testada por simulação numérica 
num programa em linguagem ASSEMBLER, executado 
num computador HP2115, do Instituto de Biologia Gul- 
benkian. 


119 


INÍCIO 


ii 


4 ==2 
+ çã | 


1d,3d | 
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Fig. 1 — Fluxograma para conversão infix-postfix (4 7) 
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Fig. 2 — Diagrama de estados de uma máquina para conversão infix-postfix (5 8) 


10. Redução dos estados da máquina [1] 


Os pares de estados compativeis-1 são (q, q), 
(o, Gu), (o Ga), (go, qu), (q q:), (q qul.las q:), (qu. qu), 
(qu Qu) (qu q), (as qo), (as. ql, (qe, qu. 

Os conjuntos máximos de estados compativeis-1 
são (Quo Qu Qu Gr), À) Qu Ga Qu), la q), (au q), tas, 
Ge Q:). Note-se que, se o par (q: q:) fosse compativel-1, 
seria possivel reduzir o número destes conjuntos. 

Resolveu-se tentar alterar a máquina de modo a que 
q e q: fossem compativeis-1, tentando manter as com- 
patibilidades entre estados já verificadas. Para que q: 
a q; sejam compatíveis-1 é necessário alterar a máquina 
de modo a que sejam também iguais as saídas corres- 
pondentes à aplicação de elementos de “na q e q. 
Uma tentativa neste sentido conduziu às tabelas Ill e IV; 
as alterações correspondem a registar o simbolo ( na 
fita auxiliar, também nos casos em que se segue a uma 
operação unária. Note-se que foi possível eliminar o 
estado Qu. 

Agora os pares de estados compativeis-l são 
(qo, Qu), (Qu Qul, (Oo Q:), (au ql, (au qul, (qu qu), (qu, 
Qu). (qu qu), (qu qu), (as, qr). 

Os conjuntos máximos de estados compatíveis-1 
são By — [Qu Qu Qu Qi), Bu — la, Qu), Bu = iq qua, 
Bu = [qu q). 

Analizando a tabela V, podemos concluir que, des- 
tes conjuntos, só B, contém estados que não são com- 
patíveis-2: q; e q. Podemos então concluir que 
os conjuntos máximos de estados  compativeis- 
-2 são Bu = (Go, Qu GQ, GR), Ba la, Q), Ba to qi), Bus 
= [qu q. É fácil verificar que estes conjuntos são tam- 
bém os conjuntos máximos de estados compativeis-3. 
Assim, a classe final da máquina é :=[Ba, Bo, Ba, 
Ba). 

Existem várias classe C mínimas, Optaremos por 
C — [Co Ci Co Ca), onde C4= iq] qua), C= 
= [qi), C.= [a], Cs = [g:), nois é conveniente que o es- 
tado inicial q, e final q; da máquina estejam no mesmo 
elemento da classe C, para que não seja necessário 
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reiniciar o estado cada vez que a máquina é operada 
de modo a traduzir outra expressão. 

A máquina classe C é então uma máquina de Tu- 
ring com 3 fitas T.=(1,8,,0,8,,% Co), onde: 


— o conjunto A, referido na definição (85), é o 
mesmo da máquina T (89), 


E S, meio 
— as funções 5, 0 À, são definidas a partir de 


duas funções õ; E x; + de um subconjunto 
dae : X C respectivamente em C e O, tais que 


FA ( 4. Q;) = 8c(x,q;) sse Xe é 
A tE g)=A tea) sse xe: 
A partir das tabelas anteriormente referidas é pos- 


sível construir as tabelas Vl e VIl, definidoras de 8 


e de à, respectivamente, 

Esta máquina também foi testada por simulação nu- 
mérica num programa em linguagem ASSEMBLER, exe- 
cutado num computador HP2115, do Instituto de Biolo- 
gia Gulbenkian. 
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Uma Bibliografia de APL 


SUMÁRIO 


Este artigo é uma colectânea de breves criticas bi- 
bliográficas de livros sobre a linguagem APL. A varie- 
dade de tópicos tratados nestas obras reflecte a pro- 
funda importância assumida por esta linguagem em 
todos os domínios da informática. 


Há mais de quatro anos, num trabalho apresentado 
no Congresso Hispano-Luso de Informática, chamâámos 
a atenção para a importância que a linguagem APL po- 
deria vir a assumir quer como notação matemática formal 
no desenvolvimento e expressão de algoritmos de cálculo 
quer na programação de computadores ou, até mesmo, 
no desenvolvimento dos próprios computadores. 

Desde então, e sobretudo a partir de fins de 1973, 
todos estes aspectos têm passado por grande expansão, 
em particular o último, agora marcado pelo aparecimento 
de micro-computadores APL comercializáveis, conhecen- 
do-se actualmente pelo menos três: um de fabrico cana- 
diano, outro americano e o terceiro franco-americano, 
não falando já no projecto do célebre STAR da CDC. 

Facilmente acessíveis em custo, estas máquinas 
podem, por outro lado, vir a ser as pioneiras de uma 
nova era da informática de gestão e administração, como 
o APL tem sido no domínio do cálculo para fins cienti- 
ficos. 

Ao mesmo tempo a bibliografia de obras sobre 
esta linguagem e as suas aplicações tem-se enriquecido 
muito. Este facto, aliado à constatação de o APL ser 
ainda, infelizmente, quase desconhecido entre nós, bem 
como de ter sido recentemente objecto de profundos es- 
tudos e aperfeiçoamentos, levou-nos à elaboração deste 
artigo, de índole estritamente bibliográfica, onde são 
apresentados alguns livros sobre APL de carácter acadé- 
mico, acompanhados cada um de uma pequena crítica. 

Esta lista não é, obviamente, exaustiva e já de si 
a escolha das obras aqui em revista reflecte um certo 
carácter do subjectividade; também o conteúdo e a ex- 
tensão de cada crítica deixam transparecer preferências 
pessoais. Por jeso mesmo lhe chamámos Uma Bibliogra- 
fia de APL, 

Em particular só apresentamos títulos em francês 
e inglês. Desconhecemos, aliás a publicação de livros 
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ABSTRACT 


This survey paper brings together a collection of 
brief reviows of books on APL. The manifold of topics 
dealt with in these works shows the primary importance 
achieved by the language in all fields of information 
processing. 


sobra este assunto em quaisquer outras línguas, com 
excepção de uma publicação dinamarquesa; desde já 
agradecemos a indicação de qualquer omissão. De qual- 
quer modo estamos confiantes no interesse deste artigo 
como fonto de referências para os informáticos portu- 
gueses. 


A linguagem APL tem suscitado profundas contro- 
vérsias e, tal como o ALGOL, é adorada até ao fanati- 
cismo por alguns e menosprezada por outros. Às críticas 
giram quase sempre em torno de dois aspectos: sintaxe 
a eficiência de execução nos computadores existentes 
actualmente. Quanto à sintaxe parece-nos nada dever 
acrescentar aqui: é, quase sempre uma questão de gosto, 
a o próprio APL permite, em geral, múltiplas variantes 
do mesmo algoritmo ou processamento. A eficiência de 
execução, é muito mais objectiva e um pouco mais drás- 
tica, em particular se compararmos o APL com lingua- 
gens como o FORTRAN em situações onde os programas 
compilados são executados repetidamente (mesmo assim 
numerosos autores referem casos onde a vantagem pende 
para o lado do APL e nós próprios apresentémos alguns). 

Embora nada impeça a existência de compiladoras 


APL le alguns de tipo incremental estão disponíveis, 
tendo mesmo um dos primeiros sistemas APL experi- 
mentais —- para os 7090 -—— sido deste tipo), a lingua- 
gem é, por natureza, adequada para uso interpretativo; 
a isso acentua-se naqueles interpretadores que aprovei- 
tam o carácter recursivo da definição de muitas das 
funções da linguagem e nos quais certas funções nati- 
vas são, elas próprias, escritas já em APL. É o caso, por 
exemplo, do APL/1100. 


Ora. a essa objecção querem vir agora responder 
também estes micro-computadores APL a que acima nos 
referimos e outros de maior porte, que certamente apa- 
recerão de futuro, os quais pelas suas características de 
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preço, universalidade e versatilidade acabarão por invadir 
também o nosso pais, 

As obras em análise são apresentadas por ordem 
alfabética e por anos de publicação. No caso de, da 
mesma obra, conhecermos duas ou mais edições, só nos 
ocuparemos da mais recente, geralmente a mais aperfei- 
coada. Indicações suplementares poderão ser obtidas 
junto do autor destas notas. 


Abrams, Philip et Lacourly, Gérard (1972). Infor- 
matique par Télephone. Langage de Programmation APL. 
Hermann et Cie., Paris. 


Trata-se de uma obra muito acessível, escrita com 
notável preocupação pedagógica por dois conhecidos 
cultores da linguagem, em particular o primeiro autor, 
participante nos projectos dos dois primeiros interpre- 
tadores operacionais de APL, há uma dezena de anos. 
Este livro não pressupõe quaisquer conhecimentos de 
programação ou de computadores, nem mesmo de ma- 
temática para além das operações fundamentais e dos 
rudimentos de álgebra. 

Contudo, muitas das funções nativas da linguagem 
não são discutidas e o mesmo sucede com quase todos 
os operadores (limitação grave mesmo para aplicações 
de complexidade intermédia), tendo a sua apresentação 
sido anunciada para um segundo volume que não deve, 
porém, vir a aparecer. De qualquer maneira, mesmo ape- 
sar das limitações de conteúdo, aconselhamos vivamente 
a sua leitura a quem deseje uma boa introdução à lin- 
quagem. 


Brooks, Frederick P. and Iverson, Kenneth E. (1968). 
Automatic Data Processing. System/360 Edition. John 
Wiley and Sons, New York. 


Uma das características importantes desta obra 
consiste na descrição do Sistema/360 e é, nor isso, de 
grande utilidade no estudo desta linha de computadores, 
bem como da sérias sucessora, o Sistema/370 (a pri- 
meira edição deste livro apresentava precisamente os 
7090; compare-se com lIverson (1962)). 

Não se trata pois, na verdade, de um livro sobre 
APL, mas antes de uma obra de informática utilizando 
a «linguagem de lIverson» (ver, de novo, abaixo Iverson 
(1962)), na descrição de algoritmos e de «hardwares». 

Notáveis são, na opinião do crítico, os capítulos 
sobre a representação interna da informação no compu- 
tador, bases de numeração, lógica simbólica, circuitos, 
aritmética de máquinas binárias e décimais, organização 
de computadores, ordenações e buscas. O mesmo sucede 
com os sistemas de exploração, em particular o 05/360, 
sendo de notar que F. P, Brooks foi um dos arquitectos 
desto sistema; por conseguinte a leitura deste capítulo 
será de grande interesse para quem se quira iniciar no 
OS via APL ou pretenda representar desta forma algorit- 
mos relacionados com o funcionamento de sistemas de 
exploração. A discussão sobre interpretadores, assem- 
blers e compiladores no capitulo 8 é do melhor que 
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alguma vez vimos sobre o assunto e recomendamo-la 
também vivamente. 

Este livro contém ainda numerosos exercícios de 
bom nivel e, ao pretender resolvê-los, mesmo o progra- 
mador experiente tirará deles numerosos ensinamentos. 
Muitos destes exercícios relacionam-se com os tópicos 
que acima passámos em revista e alguns prendem-se com 
motivos de programação e, como tal, poderão ser abor- 
dados em APL actual ou qualquer outra linguagem (no- 
tamos aqui, de passagem, que os autores repudiam a 
designação de «linguagem de alto nivel», aspecto onde 
concordamos em absoluto com eles). 


Bucklay, John W., Nagaraj, Mallur R., Sharp, Dur- 
win L. and Schenck, James W. (1974). Management 
Problem Solving with APL. Melville Publishing Company, 
Los Angeles, California. 


Um livro muito bem feito, visando o uso da lingua- 
gem para fins de gestão, com múltiplas aplicações em 
aspectos tão diversos como contabilidade empresarial, 
aplicações financeiras, estatística e métodos matemáti- 
cos na empresa, incluindo investigação operacional, ges- 
tão de existências, etc. Enfim, parece-nos que este livro 
tem, além de outros, o mérito de destruir o mito de o 
APL ser adequado unicamente para fins científicos, mito 
talvez originado pelo facto de ser uma linguagem muito 
rigorosa do ponto de vista da sintaxe e por isso exigir 
uma certa familiaridade com o instrumental matemático, 
e ter sido até agora usado quase só a partir da termi- 
nais em sistemas de tempo repartido. 

No livro, além da descrição do modo de funcion-- 
mento dos programas e exemplos da sua utilização, são 
listadas quase todas as funções que integram esses 
programas; porém, em geral, não estão optimizadas e por 
vezes constituem até maus exemolos de codificação em 
APL. Isso pode vir a tornar-se um aspecto negativo do 
livro, dada a esfera de interesses dos leitores potenciais, 
estando aind), a nosso ver e ao invés da oninião dos 
autores expressa na introdução, a leitura dependente de 
uma boa base de experiência em programação APL para 
so tirar partido das funções apresentadas, para as re- 
escrever é até completar, Neste género de aplicações um 
ajustamento a cada caso concreto é sempre impres- 
cindível. 

Esperamos que estas imperfeições, perfeitamente 
admissíveis numa obra pioneira neste vasto campo, pos- 
sam ser corrigidas em edições ulteriores. De qualquer 
modo este livro será um marco no desenvolvimento da 
linguagem, em particular se os nóveis micro-computado- 
res APL vierem a ter projecção no seu uso para outros 
fins que não os científicos. 


Demars, G., Rault, J.-C., et Ruggiu, G. (1974). Le 
Lengaçe ei les Systêmes APL. Masson et Cie., Par's. 


Talvez a única obra, de todas aqui em crítica, que 
permite uma ampla visão das possibilidades actuais do 
APL, pois cobre um largo espectro de implementações, 
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desde o clássico APL/360 e do recente APLSV (versão 
do memória virtual, muito aperfeiçoada, para os Siste- 
mas/370) da IBM, à versão APL * PLUS da Scientific 
Time Sharing Corp., e às da Burroughs, DEC-PDP, Cll, 
UNIVAC, HIS-GE e CDC. 


Trata-se de um livro de leitura por vezes bastante 
dificil, com excelentes exemplos das possibilidades mais 
avançadas da linguagem (breve introdução ao acesso a 
ficheiros, funções de execução, formato e conversão, 
etc.). Embora não pressuponha quaisquer conhecimentos 
prévios, na nossa opinião não é acessível ao principiante, 
destinando-se mais a conhecedores experientes. 


Aliás, a maneira como certas funções são apresen- 
tadas e até mesmo interligadas (lembramo-nos, por 
exemplo, da definição da reflexão via subida de compo- 
nentes ao longo da coordenada escolhida, da rotação 
via resíduo, etc.) obriga a um domínio perfeito da se- 
mântica de todas as funções da linguagem e isso só 
com longa maturidade é viável; técnica semelhante é, 
em parte, usada em Robinet (1971), como veremos 
adiante. Neste livro certas funções e operadores ainda 
não implementados mas já definidos na teoria (como 
sejam subida e descida aplicadas a quadros de número 
dimensional igual a dois ou superior) são discutidos 
sumariamente, 

No estudo dos operadores a introdução do varri- 
mento acaba por conduzir à determinação dos elementos 
neutros das funções diádicas nativas, onde para cada 
uma o elemento neutro à direita, caso exista, é precisa- 
mente o resultado da redução do vector vazio por essa 
função; se não existir, tal redução não é possível, po- 
dendo contudo estar ainda definido o varrimento para 
quadros por essa mesma função. Note-se que os ele- 
mentos neutros das funções de máximo e mínimo são, 
respectivamente, o menor e o maior valores codificáveis, 
em virgula flutuante, no computador e por isso depen- 
dentes da máquina, na realidade uma ambiguidade ainda 
sem solução no APL. 

Tal tipo de apresentação é uma raridade em traba- 
lhos deste âmbito, mais ao estilo de artigos altamente 
especializados. O mesmo sucede também quanto à recur- 
são e utilização da função de Ackermann, nesse caso 
para classificar as linguagens de acordo com a ordem 
máxima da função directamente executável (dizem os 
autores que, enquanto o ALGOL 60, o FORTRAN e o 
PL/|I não passam da ordem três, o APL permite exprimir 
a função de Ackermann de qualquer ordem). Estes e ou- 
tros tópicos do género serão, por certo, de grande inte- 
resse nos desenvolvimentos futuros. 


Quanto a estes desenvolvimentos futuros os autores 
apontam várias linhas, algumas das quais sabemos es- 
tarem já a ser atacadas com êxito por vários investiga- 
dores: introdução na linguagem de quadros generaliza- 
dos (também chamados quadros de quadros), com ex- 
tensão do domínio das funções nativas actuais e defini- 
ção de outras, estruturas de controlo mais gerais, amplia- 
ção do âmbito das funções definidas e sua extensão a 
argumentos dos operadores (em particular aos de pro- 
duto interno, produto externo, redução e varrimento) e, 
finalmente, máquinas APL, assunto a que nos referimos 
logo no início deste artigo e esteve na sua origem. Vol- 
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taremos a abordar todos estes aspectos das tendências 
actuais da investigação no domínio do APL quando, 
abaixo, fizermos, a análise sumária das actas dos con- 
gressos realizados na Europa em 1971 (Paris), 1973 
(Copenhaga) e 1975 (Pisa). 

Resumindo: para o leitor com conhecimentos prá- 
vios a prática da linguagem esta obra é essencial para 
fazer comparações das várias versões referidas acima 
o efectuar trabalhos nos domínios de investigação se- 
guidos actualmente. O livro conclui com uma extensa bi- 
bliografia de dezenas de títulos, livros e artigos de re- 
vistas, englobando numerosos trabalhos de investigação 
nesta campo, e tornando assim a obra ainda mais va- 
liosa. 


Dupont, Philippe (1973). Mathématiques pour Infor- 
maticiens de Géstion. Société d'Enseignement Supérieure, 
Paris. 


Este não é, de facto, um livro de APL, e o titulo 
traduz bem o seu conteúdo. Porém, um dos seus capitu- 
los é dedicado à «linguagem de lIverson» (veja-se Brooks 
and Iverson (1968) e lIverson (1962)), ainda sob a 
forma definida originariamente e, por esta ter sido a an- 
tecessora do APL, o incluímos aqui. 

Esse capitulo integra-se, aliás, no conteúdo e na 
orientação do livro. Como todos os outros capítulos é de 
leitura muito agradável, cobrindo o conjunto uma boa 
dose de conhecimentos de matemática necessários à 
bagagem de qualquer informático, mesmo para os que 
se ocupam com problemas estritos de direcção e aplica- 
ções administrativas, 


Dupont, Philippe et Tallineau, Yves (1974), Les 
Langages de Programmation en Parallele: FORTRAN, CO- 
BOL, PL/l et APL. Masson et Cie., Paris. 


O titulo é bem explícito do conteúdo do livro: uma 
colectânea de bons problemas programados nas quatro 
linguagens de maior interesse prático nos nossos dias. 
Dizem os autores, com razão, que um informático conhe- 
cedor de uma ou duas linguagens de programação, dese- 
joso de compreender a lógica, os objectivos e um pouco 
da semântica de uma outra se encontra invariavelmente 
na situação de um adulto quando pretende aprender uma 
nova língua. Também aqui a tendência é para os novos 
conhecimentos serem adquiridos por recurso e referência 
aos anteriores. 

É esta a linha onde se enquadra o trabalho de Du- 
pont e Tallineau. Nos exemplos aí discutidos o APL apre- 
senta, quase sempre, vantagens nitidas. Os seis casos 
apresentados vão de simples cálculos e processamentos 
do matrizes e quadros a chamadas de subprogramas e 
aplicações não-numéricas, cobrindo um vasto domínio 
de possiveis utilizações de cada uma das linguagens. 

Como aspecto negativo no que se refere ao APL, 
as funções desta linguagem estão codificadas em APL/ 
/1130, precisamente a primeira versão operacional da 
linguagem (divulgada em 1968) e por isso muito ultra- 
passada hoja em extensão, flexibilidade e rapidez de 
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execução, devendo o leitor interessado ter em considera- 
ção as variantes profundas existentes agora, Aliás, O 
mesmo sucede no caso do PL/I, para o qual os autores, 
sem qualquer explicação, usaram também um compila- 
dor antigo. 

Em edições ulteriores deste livro, o conteúdo deve- 
ria ser por este motivo bastante actualizado; uma remo- 
delação tornaria a obra de utilidade ainda maior, quer 
para o informático orientado para questões práticas ou 
até mesmo para o gestor, interessados sobretudo em li- 
rar partido concreto do que individualmente elas permi- 
tem realizar. 

Uma boa parte deste volume é dedicada ainda a um 
breve estudo da sintaxe das linguagens e a técnicas de 
programação especificas de cada uma delas. Somos 
também da opinião ser este capitulo susceptivel de me- 


lhoria, aumentando a utilidade global do livro, e espera- 
mos que isso se venha a verificar no futuro. 


Falkoff, Adin and Iverson, Kenneth E, (1972). APL/ 
/360 OS and APL/360 DOS Users Manual. Internatio- 
nal Business Machines Corporation, White Plains, New 
York, Form No. H20-0906. 


Esta é a verdadeira «bíblia» da linguagem, oficial- 
mente a definição correcta e exaustiva do APL/360, por 
assim dizer o «relatório revisto do APL» (o livro de Pa- 
kin (1972) sendo, nesta terminologia, o «relatório não 
revistos). É constituido por uma exposição muito clara, 
não obstante conter alguns exemplos muito complexos, 
compactos e difíceis de acompanhar. 

Como definição oficial do APL/360 engloba uma 
apresentação completa de todas as funções e operadores 
da linguagem e é por isso indispensável para quem tenha 
do usar este sistema, sendo de notar terem as várias 
edições (esta é a terceira e, a avaliar pelo anúncio em 
1973 do APLSV, talvez a última) saido sucessivamente 
aumentadas, não obstante serem decalcadas da primeira, 
de 1968. 

Os apêndices B e € são, sem dúvida, os de maior 
profundidade, neles se tratando de exemplos avançados 
o de novas funções da linguagem, até à última versão 
disponivel do APL/360 (por vezes designada vor XM6, 
do código do programa-produto). Este manual abrange 
ainda um tratamento detalhado de todas as operações 
de terminais e dos comandos do sistema, embora nem 
uns nem outros sejam inerentes à linguagem. 

Esta publicação é absolutamente essencial, até por 
ser wc manual do construtor, para quem utilizar o APL/ 
/360. Por isso ficaremos por estes comentários resumi- 
dos, guardando outras notas para as críticas dos restan- 
tes livros de lIverson. 


Falkoff, Adin and Iverson, Kenneth E. (1973). The 
APLSV User's Manual. International Business Machines 
Corporation, White Plains, New York, Form No. H20- 
-1460. 


Trata-se, por assim dizer, do complemento da obra 
acima, estendendo-se a definição da linguagem ao APL- 
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SV (APL Shared Variables) a que fizémos já várias re- 
ferências) e, na generalidade são ainda aqui válidas as 
mesmas considerações. 

É de notar que nem todas as funções trabalham 
exactamente da mesma maneira nas duas versões. Em 
particular, fornecem resultados distintos as de resíduo, 
codificação e descodificação, todas diádicas, e a trans- 
posta monádica de quadros de número dimensional su- 
perior a dois. Esta última generaliza agora a quaisquer 
quadros o conhecido teorema da álgebra linear de ser 
a transposta do produto de matrizes o produto por or- 
dem inversa das transpostas dos factores, verificadas as 
condições de conformabilidade. 

Esta versão do APL é muito mais completa do que 
quaisquer outras anteriores, nela estando englobadas to- 
das as novas possibilidades da linguagem já definidas 
do ponto de vista teórico (com raríssimas excepções): 
acesso a ficheiros, «interface» com outras opções do sis- 
tema de exploração (05/VS) e compiladores de todas 
as outras linguagens, diálogo (via sistema de exploração) 
com todas as unidades periféricas, variáveis partilhadas 
(as Shared Variables), novas funções de formato, con- 
versão o execução e o operador de varrimento, etc. Este 
manual é, de parceria com o acima, de consulta fre- 
quente e obrigatória para todos os que tiverem de usar 
o APLSV, e ainda para todas as pessoas que queiram 
ler a literatura recente sobre a linguagem e suas aplica- 
ções ou enveredar pela pesquisa. 


Freiberger, Walter and Grenander, UIf (1971). 4 
Course in Computational Probability and Statistics. Sprin- 
ger Ve lag, New York-Heidelbe g. 


Este é um livro de estatística e probabilidades (para 
sermos mais concretos diriamos, de probabilidades e 
processos estocásticos), com bons exemolos de aplica- 
ções práticas em APL. Não é, concerteza, acessivel a 
quem não domine o instrumental probabilístico ao nível 
de um curso intermédio, e é por isso de menor interesse, 
por exemplo, para engenheiros e economistas, mas de- 
verá sem dúvida ser de agrado para matemáticos e fi- 
sicos. 

Algumas das aplicações do APL aqui descritas são 
muito interessantes e pouco vulgares: processos de nas- 
cimento e morte, crescimento de populações, processos 
de Markov, simulação, etc. Pena é os requisitos mate- 
máticos serem elevados e por esse facto o livro não ter 
tido a divulgação que merece e o prestígio dos seus 
autores avalisa, 


Gilman, Leonard and Rose, Allen J. (1974). APL; 
An Interactive Approach. Second Edition. John Wiley 
and Sons, New York. 


Na opinião do crítico uma das melhores obra sobre 
APL para quem nada sabe e deseja possuir um bom do- 
mínio da linguagem, especialmente indicado para auto- 
estudo, facilitado se existir disponivel um terminal. Tra- 
ta-se da segunda edição de um livro (a nosso ver a pri- 
meira era já de si excelente), aperfeiçoada de modo a 
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trater com detalhe os sistemas APLSV e APL * PLUS. 
Evita dificuldades de natureza matemática, vamos lá, que 
ultrapassem os conhecimentos do 1.º ano das nossas 
universidades; mas mesmo assim chega a discutir tópi- 
cos como processamentos de quadros de várias dimen- 
sões e números dimensionais, funções de argumentos 
booleanos, inversão de matrizes e algumas variantes e 
aplicações, como, por exemplo, ajustamentos pelo mé- 
todo dos menores quadrados, spinores, etc. 

Também muitas das possibilidades da linguagem 
mais específicas são também apresentadas com detalhe, 
em particular as relacionadas com processamentos não- 
-numéricos, incluindo consultas a ficheiros e <ua geração, 
bem como a opções gráficas em terminais com impressão 
(não para as pantalhas de raios catódicos). 

Ão longo de numerosos exercícios são ainda desen- 
volvidas múltiplas aplicações, nem todas simples, e nal- 
guns casos o mesmo problema é resolvido em vários 
capítulos, do maneira a exemoplificar as possibilidades 
das várias funções nativas da linguagem e de outras 
funções definidas no texto ou deixadas para o leitor 
definir. 

De grande interesse é o capítulo sobre a codifica- 
ção interna dos valores numéricos e literais nos Siste- 
mas/370 e a explicação de certas contradições aparen- 
tes de alguns resultados (como seja a comparação lógica 
de dois valores distintos fornecer como resultado 1 se 


a distinção fôr provocada por bits de ordem inferior, o 
que, no entanto pode ser evitado por comando ou fun- 
ção do sistema). E, por termos tocado neste ponto, indi- 


camos esto livro em especial naqueles casos em que se 
queira tirar partido da possibilidade de substituir os 
comandos do sistema (estranhos à linguagem, como 
frizámos em Falkoff and Iverson (1972)) pelas novas 
variáveis e funções de sistema — muito mais flexíveis 
a susceptíveis de ser incorporadas em funções definidas, 
devendo prevenir que algumas dessas funções de sis- 
tema são exclusivas do APL * PLUS. 


Grey, Louis D. (1973). 4 Course in APL/360 with 
Applications. Addison-Weslay Publishing Company, Rea- 
ding, Massachusetts. 


Eis um livro destinado a quem pretender aplicar a 
linguagem tendo em vista questões de análise numé- 
rica, a nosso ver altamente indicado vara engenheiros. 
Apresenta óptimos problemas de matemática (mesmo 
da questões mais avançadas, como equações com deriva- 
das parciais) e estatística; tem como pré-requisitos para 
a compreensão destes exemplos bastantes conhecimentos 
desses dominios. 

Também certos aspectos relacionados com gráficos 
a formatos de saída são discutidos com detalhe. Como 
seria de esperar as funções apresentadas nestes casos 
são de grande complexidade e muito compactas, mas 
isso é bem uma característica das aplicações avançadas, 
por vezes levada a um extremismo condenável por certos 
autores, o que acaba por ser aproveitado pelos opositores 
do APL... Contudo o interesse que certas funções desse 
tipo apresentam é bem um justificativo como é evidente, 
para a sua inclusão aqui. Talvez seja isso que falta no 
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livro de Buckldy et al., mas, como é evidente, os âmbi- 
tos e as audiências de cada uma destas obras são com- 
pletamente distintos. 

Hellerman, Herbert (1973). Digital Computer Sys- 
tem Principles. Second Edition, McGraw-Hill Book Com- 
pany, New York, 


É um livro sobre computadores baseado em APL, e 
acima de tudo, um bom exemplo das possibilidades da 
linguagem na definição de sistemas. Para além de uma 
óptima, embora curta, introdução ao APL, apresenta tó- 
picos de profundo interesse, sobre o ahardware» e o 
«software» de várias máquinas, com exemplos. Assim, 
lá aparecem assuntos tão diversos e importantes como 
representação interna de dados, estudo de instruções- 
-máquina e linguagem Assembler, vírgula flutuante, es- 
crita de compiladores e modelos de tradutores, ordena- 
ções, etc. Uma série impressionante de informações, por 
certo só viável devido ao uso sistemático do APL. 

Também são apresentados vários exemplos relacio- 
nados com a construção dos circuitos internos dos com- 
putadores e a sua lógica, bem como de unidades perifé- 
ricas. Um capítulo á dedicado à arquitectura de compu- 
tadores e micro-programação, outro a máquinas de Tu- 
ring. Qualquer dos dois, limitado em extensão dado o 
carácter geral do livro, constituem magníficos exemplos 
desta via de ataque a problemas concretos. 

Para além do interesse deste trabalho, derivado do 
uso extensivo do APL, consideramo-lo como um dos me- 
lhores de carácter introdutório sobre computadores. La- 
mentamos apenas não poder, dado o seu âmbito mais 
geral o a natureza sumária das críticas desta bibliografia, 
desenvolver aqui uma análise extensiva, como seria 
nosso desejo. 


Hurtubise, Rolland et Poulin, Yves (1973). APL et 


les Cadres. Les Presses de [Université du Québec, Qué- 
bec. 


Um outro livro, talvez mais orientado para fins 
simples de gestão, elementar quer nos requisitos mate- 
máticos, quer no âmbito da linguagem desenvolvida. 
Todos os exemplos são também muito acessíveis, apre- 
sentando contudo um apêndice sobre bancos de dados, 
retratando aliás uma aplicação concreta desenvolvida 
com a linguagem na Universidade do Québec para ges- 
tão estudantil. Os bancos e bases da dados serão con- 
certeza uma outra área de desenvolvimento, ainda aqui 
a dar os primeiros passos, o que do ponto de vista his- 
tórico se compreende por, até há bem pouco tempo, o 
APL ter sido bastante limitado nas entradas e saidas e 
no acesso a ficheiros. 

Em termos de comparação colocariamos este livro 
ao nivel de Abrams et Lacourly (1972), examinado acima 


Iverson, Kenneth E. (1962). A Programming Lan- 
guage. John Wiley and Sons, New York. 


O livro do «pair do APL, o marco oficial do nasci- 
mento da linguagem (se bem que a sua génese possa ser 
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traçada a alguns anos antes), além de ter sido o seu 
título o originador do acrónimo da linguagem, por vezes 
denominada nos primeiros anos após a publicação desta 


obra, «linguagem de Iverson» (tal denominação caiu hoje 
em desuso). 


Em si a notação algébrica apresentada neste livro 
tem já muito poucas afinidades com o APL actual, mas 
muitas das ideias básicas das funções nativas definidas 
nas versões actuais e numerosos algoritmos encontram- 
-so al, pela primeira vez, com rigor notável. 

Um capítulo muito interessante deste livro descreve 
a arquitectura dos antigos sistemas 7090 e faz lembrar 
um artigo publicado em fins de 1964 por Falkoff, Iver- 
son e Sussenguth sobre a arquitectura do Sistema/360, 
hoje considerado um clássico sobre o assunto, utilizando 
também o APL na simbologia usada neste livro. Outros 
capítulos tratam de numerosos tópicos de informática, 
sendo a nosso ver particularmente interessantes os que 
se referem a métodos de ordenação (um tónico de parti- 
cular favoritismo nestes livros), árvores e grafos, subro- 
tinas, metaprogramas e microprogramação, ainda hoje 
de leitura muito proveitosa. O livro tem um pouco o 
mesmo sabor de Brooks and lIverson (1968), mas mais 
avançado, embora não tão actualizado. 


Contudo, somos dos que pensam ser o grande mé- 
rito deste livro, pelas inesgotáveis consequências que 
têm tido, a introdução de funções de natureza escalar es- 
tendiveis a quadros, elemento por elemento, e sobretudo 
a definição de alguns dos actuais operadores do APL, a 
saber, os de produto interno, produto externo, redução 
a varrimento. Com excepção do último já todos faziam 
parte das versões anteriores ao APLSV, embora ainda 
agora não estejam completamente desenvolvidas as ba- 
ses teóricas do operador de varrimento e as suas implica- 
ções em questões práticas, sobretudo no caso de se epli- 
car a funções definidas pelo utilizador. 


Como dissémos antes, muito pouca semelhança 
existe entre a simbologia matemática desenvolvida neste 
livro o o APL actual, Contudo, como introdução teórica 
ao estudo da linguagem é de importância fundamental, 
nele se vislumbrando muitas das linhas de investigação 
seguidas hoje. 


Iverson, Kenneth E. (1966). Elementary Functions. 
Scienco Research Associates, Chicago, Illinois, 


No prosseguimento do seu livro de 1962 lIverson 
escreveu esto compêndio de matemática ao nível aproxi- 
mado do nosso curso complementar dos liceus, onde, 
sem perder de vista os fins computacionais,discute na 
orientação da filosofia do APL as funções mais importan- 
tea da matemática pré-universitária: trigonométricas, ex- 
ponencial, logarítmica, polinomiais, inversas, hiperbóli- 
cas, etc, O conceito de produto interno, de importância 
basilar na linguagem, como acima vimos, é introduzido 
a partir da vulgar função polinomial (soma de produtos) 
o depois generalizado. É também muito bom e de recorte 
original, mesmo nos dias de hoje, o capítulo dedicado 
ao estudo das derivadas. Aliás pensa-se que a linguagem 
ainda deverá vir a abarcar um operador especial de deri- 
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vação e outro de primitivação, tópicos seguidos por um 
ou outro investigador, sem grandes resultados visíveis, 
de momento, 

Será de interesse mencionar aqui de passagem, 
pelo facto de ser uma característica marcante da lingua- 
gem, terem sido neste livro apresentadas pela primeira 
vez algumas justificações da célebre ordem de execução 
das operações, onde para todas as funções, sem qualquer 
excepção, é sempre da direita para a esquerda (por isso, 
a título jocoso alguns opositores chamam ao APL «ara- 
bic programming language» — outra anedota admite a 
hipáóteso de Iverson ser canhoto...). Esta característica 
só se estabeleceu quando estavam a ser escritas as pri- 
meiras implementações e é radicalmente diferente da 
ordem de execução do livro de 1962, onde ainda eram 
utilizadas as regras da álgebra. Ora um dos exemplos 
apresentados por Iverson para justificação da regra «da 
direita para a esquerda» é o cálculo do valor de fracções 
continuas regulares e de polinómios. Aliás o número e 
é mesmo introduzido via aproximação de fracção con- 
tinua, logo um exercício dum dos primeiros capítulos; 
estes a outros exemplos demonstram a profundidade 
do livro, não se devendo esquecer estar dedicado a um 
nível abaixo do universitário. 

E, como mencionámos a ordem de execução das 
funções, gostariamos de salientar ser agora quase dez 
anos após o seu estabelecimento, uma das orientações 
actuais antevistas para a linguagem a extensão desta re- 
gra a todos os operadores, considerando-se em APL como 
operador toda a função (monádica ou diádica) que pos- 
sui uma ou duas outras funções por argumentos (que 
por sua vez terão ainda argumentos numéricos). Nalguns 
casos, especialmente no da indexação para os operado- 
res de eixos, já referidos ao analisarmos ecima o livro 
de Demars, Rault et Ruggiu (1974), isso tem levado a 
uma redefinição profunda das estruturas dos quadros; 
compare-so com a breve resenha ao artigo de Lathwell 
e Mezei sobre a descrição formal do APL no exame, 
abaixo, das actas do colóquio de 1971 em Paris. 

Mas, voltando ao livro de Iverson de 1966: quando 
foi publicado existia já uma versão experimental do APL/ 
[360 e a ela é dedicado o último capítulo do livro. Con- 
tudo esta versão é anterior ao grande incremento do APL 
como linguagem de programação e, aqui o além, o âm- 
bito dessc capítulo está um pouco ultrapassado. Hoje 
esto livro seria mais um compêndio de matemática, como 
dissémos, e assim é perfeiiamente actual e mesmo 
original, mais do que propriamente um livro de APL no 
sentido moderno, ou então uma óptima introdução à 
linguagem para o ensino liceal. Contudo, um dos seus 
méritos é ainda ser de grande relevo para a compreensão 
da evolução da linguagem, em especial durante os pri- 
meiros anos da sua existência, 


Iverson, Kenneth E. (1972). Algebra: An Algorith- 
mic Treatment. Addison-Wesley Publishing Company, 
Reading, Massachusetts, 


No domínio a que se dedica, ou seja, o ensino da 
álgebra elementar com a utilização de técnicas compu- 
tacionais e do APL, este livro é, na nossa opinião, ex- 
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traordinário e um ponto de referência no campo da lin- 
guagem e do seu uso como veículo de comunicação de 
algoritmos. 

São tratados tópicos tais como a resolução de sis- 
temas do equações lineares, funções lineares, polinómios 
o operações, avançando mesmo por caminhos desconhe- 
cidos ao nivel da nossa matemática liceal, como sejam 
operações elementares com matrizes e uso da inversa 
para a resolução de sistemas de equações, ajustamentos 
pelo método dos menores quadrados (via função nativa 
de inversão de matrizes, mas sem apresentar o prócesso 
de cálculo, como seria de esperar); o mesmo sucede, 
eliás com o conceito de determinante, definido, sob 
certas condições, como função de matriz quadrada. 

A definição do produto externo é introduzida muito 
cedo, logo no início do capítulo 2, em conexão com a 
construção de tabelas das funções soma e produto. 
O produto interno aparece também bastante mais cedo 
do que seria de -sperar, e é curiosa a sua introdução 
num algoritmo de caminho crítico, muito simples como 
conviria para a audiência do livro. A redução surge de 
forma ainda mais natural, mesmo no princípio do livro, 
mal está definida a soma de dois números inteiros, na 
página 7. 

Um capítulo relaciona vários aspectos de aritmética 
p representação dos números, em especial dos números 
racionais; é dos mais curiosos e merece estudo detalha- 
do. Os números primos tinham sido antes apresentados 
por via programática e aí usados na definição do menor 
múliiplo comum e do máximo divisor comum. Outros 
assuntos tratados de maneira notável são os coeficientes 
binomiais, polinômios factoriais e indução matemática. 

Tanto o aprendiz de APL como qualquer pessoa in- 
teressada nas modernas exposições de matemática de- 
viam considerar este livro como uma valiosa leitura obri- 
gatória. 


Katzan, Harry (1970). APL Pregramming end Com- 
putos Techniques. Van Nostrand-Reinhold, New York. 


O título é um tanto enganador, pois além de três 
capítulos (elementares) sobre APL, trata também de mui- 
tos outros assuntos, sem qualquer ligação na exposição, 
como sejam, generalidades sobre computadores, siste- 
mas de exploração (incluindo memória virtual) e lingua- 
gem FORTRAN. Está na mesma linha de outros livros 
desto autor e deve ser lido nesse contexto. 

A nosso ver poderia ser usado ao nível de cursos 
de introdução aos computadores das nossas universida- 
des, dado o seu carácter de obra de referência. O apên- 
dice é também dedicado ao APL e, aí sim, apresenta 
bons exemplos, alguns respigados da estatística, investi- 
gação operacional e análise numérica, O que aumenta 
substancialmente a sua utilidade. 


Katzan, Harry (1971). APL Users Guide. John 
Wiley and Sons, New York. 


Esto manual é, em boa parte, uma reposição dos 
capítulos sobre APL (com excepção do apêndice) do 
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livro acima analisado, com exemplos talvez melhorados. 
Pode também servir como obra de introdução à lingua- 
gem, mas quanto a nós a sua grande utilidade reside na 
extensa análise dos comandos do APL/360, operação 
dos terminais da linguagem e mensagens do sistema, 
em especial das mensagens de erro. 

Quando estes dois livros de Katzan foram publica- 
dos havia relativamente poucas obras sobre estes assun- 
tos de carácter didático. Desempenharam, por isso, um 
papel importante, embora hoje estejam um pouco ultra- 
passados. 


von Maydell, U. E. (1970). An Introduction to Pro- 
bability Using APL. Department of Computing Science, 
Tho University of Alberta, Edmonton, Alberta. 


Embora seja uma publicação de características um 
pouco especiais, como o seu título indica, pensamos 
mesmo assim que deve figurar numa colectânea de cri- 
ticas deste género, especialmente como prova da versa- 
tilidade da linguagem. 

Nesta manual a autora preocupa-se com o ensino 
das probabilidades a um nivel introdutório, e isso é feito 
via APL. Por exemplo, o conceito de valor esperado de 
uma variável aleatória discreta assumindo os valores 
quo constituem o vector X com probabilidades P pode 
ser escrito em APL muito simplesmente +/PXX ou, 
como uma alternativa para o uso do produto interno, 
P+H.XX. 

Esta técnica é hoje seguida por numerosos autores, 
mesme no campo da estatística e das probabilidades. 
Pensamos nos livros abaixo de Smillie (1973, 1974), 
numas notas, já há anos aguardando publicação, de Ans- 
combe (Statistical Computing with APL), e noutros tra- 
balhos, mormente artigos de revistas. O texto de Frei- 
berger and Grenander (1971) vai numa orientação um 
pouco diferente. 

Contudo o valor do trabalho de von Maydell reside 
no facto de ter sido um dos trabalhos pioneiros. A publi- 
cação está bastante limitada, tendo sido anunciado um 
segundo volume, infelizmente ainda por publicar, pois 
pensamos que trabalhos deste tipo podem desempenhar 
um papel de relevo especial na projecção da linguagem. 


Mock, Theodore J. and Vasarhely, Miklos A. (1972) 


APL for Management. Wiley-Becker and Hayes, New 
York. 


O titulo do livro é bem representativo do seu con- 
teúdo, tratando temas de gestão e planeamento, um pou- 
co na linha de Buckley et al. (1974). Discute tópicos 
tais como estatística dedutiva elementar (vejam-se os 
livros de Smillia abaixo), programação matemática, simu- 
lação, PERT. etc. além de uma apresentação simplifi- 
cada das funções nativas mais usadas na linguagem, sem 
se preocupar com processamentos de quadros de número 
dimensional para além de dois; tendo em vista as fina- 
lidades do livro tal restricção não é limitação alguma. 

Muitas das funções usadas, por exemplo as da 
nrogramação linear interactiva, devem-se ao Professor 
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Smillia (cujos livros são examinados no seguimento) e 
integram-se num conjunto de programas de aplicações 
estatísticas e matemáticas conhecido por STATPACK2, 
disponível na Universidade de Alberta. 

Pensamos ser este o livro ideal para o gestor inte- 
ressado em extrair algum proveito da linguagem sem se 
preocupar com detalhes de escrita de programas, manu- 
tenção do sistema, etc. Nesta óptica recomendamo-lo 
sem hesitações, 


Pakin, Sandra (1972). APL/360 Reference Manual. 
Second Edition. Science Research Associates, Chicago, 
Hlinois. 


A primeira edição deste livro data de 1968, con- 
temporânea das duas primeiras versões disponíveis da 
linguagem para uso geral e foi por isso um dos primeiros 
que, historicamente, contribuiram para a sua expansão. 
Essa edição é anterior mesmo ao User's Manual, Fal- 
koff and Iverson (1972), referido cima e que então cir- 
culava sob forma policopiada. Por isso mesmo chamá- 
mos a este «o relatório da linguagem APLs», 

Pelos motivos acima e pela preocupação de trata- 
mento exaustivo de sintaxe da linguagem são-lhe feitas 
frequentes referências na literatura quanto à definição 
do APL. Por ter em consideração os conceitos introdu- 
zidos entre 1968 e 1971, esta segunda edição do livro 
é mais completa do que a primeira, sem contudo intro- 
duzir modificações profundas, e conservando no essen- 
cial o mesmo sabor, sempre com grandes cuidados pe- 
dagógicos. 

À discussão de todos os elementos da linguagem e 
funções nativas é muito rigorosa. Também as funções 
definidas e as técnicas da sua escrita, edição e altera- 
ção são tratadas com toda a atenção e do livro fazem 
ainda parte numerosos exemplos, em geral bem exempli- 
ficativos das enormes possibilidades da linguagem. 

Não obtante não considerar ainda a maior flexibili- 
dado oferecida pelo APLSV trata-se de um excelente li- 
vro de consulta, muito acessível em particular para quem 
tiver de utilizar na prática a linguagem através de um 
terminal. 


Penfield, Paul (1971). MARTHA User's Manual. The 
MIT Press, Massachusetts Institute of Technology, Bos- 
ton, Massachusetts. 


Um manual sobre a utilização de um conjunto de 
programas, designado por MARTHA, destinados a ana- 
lisar circuitos eléctricos de tipo linear e alguns de tipo 
não-linear. 

Não nos podendo pronunciar sobre o conteúdo 
técnico deste documento, notamos apenas que inclui 
um sumário muito ligeiro da linguagem. Estamos em 
crer acabarem muitos utilizadores destes programas por 
sentir a necessidade de aprofundar essas bases de APL, 
tendência louvável neste tipo de aplicações em que o 
computador é um precioso auxiliar de ensino e ainda 
uma fonte de simulação de experiências, Numerosos 
casos de saidas gráficas. 
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Perlis, Alan J. (1975). Introduction to Computer 
Science. Harper & Row, New York. 


O Professor Perlis não necessita de qualquer apre- 
sentação, e o seu prestígio é uma garantia da qualidade 
do livro, cujo conteúdo transparece bem no título. Trata- 
-se de outra obra de valor incalculável, por um lado pelo 
facto de ser usado não só o APL como também o ALGOL 
na descrição de algoritmos, arquitectura de computado- 
res, processamento de linguagens, organização de dados 
a análise de sistemas, Parece, pois, não serem as duas 
linguagens tão incompatíveis como alguns académicos 
pretendem fazer crer, pelo menos para este tipo de ex- 
posição e para uma autoridade como é o autor. 


Por outro lado, parece-nos interessante mencionar, 
o Professor Perlins é, de facto, um dos implementadores 
do APL, tendo sido responsável pela escrita do interpre- 
tador disponível para o DEC PDP-10 desenvolvido na 
Carnegie-Mellon University há cerca de dois anos e, 
segundo se consta, um dos melhores de qualquer sistema 
Pensamos que tal só é válido na presente lista de re- 
ferências para quatro outros autores, o Dr. Philip Abrams 
(acima mencionado em Abrams et Lacourly (1972)), o 
Dr. Bernard Robinet (ver Robinet (1971)) e Richard 
Lathwell a Jorgo Mezei (de quem um importantíssimo 
trabalho é referido no Coloque APL (1971) examinado 
abaixo). 

Fizémos, no início deste artigo, uma menção ao 
carácter interpretativo por excelência do APL. Ora, Perlis 
chegou ao ponto de em 1971 o ter já considerado como 
uma extensão natural do Assembler. É curioso como 
esta afirmação profética está agora em vias de tomar 
forma, O livro reflecte estas facetas das contribuições de 
Alan Perlis para o progresso do APL. 


Polivka, Raymond and Pakin, Sandra (1975). APL 
The Language and lts L'sage. Prentice Hall, Englewood 
Cliffs, New Jersey. 


Trata-se de um livro muito ao estilo do de Pakin 
(1972), com a colaboração de um dos primeiros cul- 
tores da linguagem, Polivka, autor de um sistema de 
ensino de APL com o auxílio do computador, escrito na 
própria linguagem e que fazia parte das primeiras versões 
do APL/360; um exemplo do uso do APL no CAI (Com- 
puter ÁAssisted Instruction). 


Prager, William (1971). An Introduction to APL. 
Allyn and Bacon, Boston, Massachusetts. 


Este é um outro compêndio muito interessante, ao 
nivel de complexidade dos de Abrams et Lacourly (1972) 
de Pakin (1972) e de Hurtubise et Poulin (1973), acima 
referenciados, talvez de complexidade intermédia entre 
os dois primeiros. Como Abrams e Lacourly também 
não apresenta todas as possibilidades da linguagem, po- 
dendo por isso servir de texto para uma introdução ao 
nível da transicção liceu /universidade ou para quem te- 
nha em vista avaliar das possibilidades concretas do APL 
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e. em casos concretos, fazer comparações com outras 
linguagens. 

Alguns tópicos tratados exigem certos conhecimen- 
tos matemáticos como cálculo integral, matrizes e deter- 
minantes, mas nunca muito profundos. Muito bem de- 
senvolvidas aqui são as saidas de gráficos, que foram 
omitidas de todas aquelas obras; pela facilidade do seu 
manuseamento em APL são de exploração assaz com- 
pensadora a este nível. 


Robinet, Bernard (1971). Le Langage APL. Editions 
Technip, Paris. 


Trata-se de um livro impecável do ponto de vista 
científico, esgotando completamente as possibilidades 
do APL/360. Enquadra óptimos exemplos, alguns exi- 
gindo conhecimentos de matemática um pouco avança- 
dos, mas para leitores com uma tal bagagem será dos 
que mais vivamente aconselhamos. O autor, professor 
no Instituto da Programação da Universidade de Paris, 
foi também um dos escritores do compilador incremental 
APL do ICL-Elliott 4130, devendo-se-lhe ainda vários 
estudos teóricos sobre a semântica dos operadores de 
quadros. 


O livro abre com um prefácio do Professor J. Ársac, 
logo um inimigo acérrimo do APL (mera coincidência 
ou atitude deliberada de Robinet?); conhecido assim o 
pensamento de Arsac, a leitura deste prefácio é uma 
imposição, quase diriamos, imprescindível, nara quem se 
dedicar à linguagem. De facto, além de um ataque ao 
contexto e à forma do APL («acabam-se os belos ciclos 
de cálculo e as formas recursivas»), Ársac antevê as 
consequências dos programas compactos («sim, a lingua 
APL é assim: diz muito com poucas palavras»), a deixa 
pairar sem resposta a questão da vantagem da expressão 
algorítmica. 


Mas, esta crítica dirige-se ao livro de Robinet e não 
ao seu prefácio. Por isso passamos adiante. Nas duas 
partes que o compõem a linguagem é dissecada, passe 
o termo, em todos os pormenores, constituindo um tra- 
tamento exaustivo do APL/360 tal como à data era de- 
finido. Por vezes a ordem porque os assuntos são trata- 
dos é um pouco estranha, embora faça sentido. Por 
exemplo, os quadros são introduzidos logo no início, 
bem como as funções que a eles se aplicam especifica- 
mente, enquanto, pelo contrário, a apresentação formal 
da ordem de execução das operações só aparece pró- 
ximo do fim da segunda parte, surgindo então como a 
única natural, até porque já foi sendo discretamente 
usada, embora só sejam apresentadas as justificações 
do próprio Iverson (1966). 


As funções definidas são discutidas em dois ca- 
pítulos distintos, um em cada parte. No primeiro é apre- 
sentado o essencial para o seu uso extensivo no decurso 
do livro, adiando-se todos os pormenores da sua escrita 
até ao outro capítulo. Notámos que também a menção 
dos comentários vó vem quase no final, opção aceitável 
embora a maioria dos autores coloque a sua explicação 
logo na primeira página do livro que escrevem, ou 
quase. 
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Como dissémos, o nível desta obra é avançado, re- 
flactindo-se nos exemplos. Para citar uma amostra: ope- 
rações sobre estruturas algébricas, processamentos não- 
-numéricos, gerações de permutações, operações sobre 
permutações e quadrados latinos ortogonais, teoria dos 
grafos, etc. 

Em conclusão, diremos que se trata de um livro 
excepcional, a par com Gilman and Rose (1974), quiçá 
um excelente aperitivo para o de Demars, Rault et Ru- 
ggiu (1974). Só é de lamentar que esteja um pouco de- 
sactualizado em relação aos desenvolvimentos da lingua- 
gem posteriores a 1971, em particular ao APLSV, o que 
esperamos venha a ser ultrapassado numa próxima edi- 
ção. 


Smiíillio, Keith W. (1973). A Short Introduction to 
APL/360. Department of Computing Sciences, The Uni- 
versity of Alberta, Edmonton, Alberta. 


Um livrinho com menos de 20 páginas de texto, 
um apêndice de pouco mais de 30 com listagens de 
exemplos e outro sobre vários aspectos elementares da 
estatística e probabilidades com menos de 10, usado 
em várias universidades do Canadá e dos Estados Uni- 
dos para introdução ao APL, dedicado a estudantes cuja 
orientação principal não é a informática, de preferência 
com alguns conhecimentos de estatística elementar. 

À selecção dos tópicos foi muito cuidada e os 
exemplos são de compreensão fácil e de grande utili- 
dade; o nivel matemático não é elevado, mas isso de 
modo algum diminui o valor do trabalho. O autor teve 
a preocupação de eliminar alguns tópicos de maior com- 
plexidade, mas os apresentados são bastante detalhados, 
mesmo a função de inversão de matrizes, quer usada 
como tal quer em ajustamentos lineares, assuntos nem 
sempre devidamente considerados a este nivel mais 
elementar. 

O apêndices 2, com as saídas gráficas usando as 
funções desenvolvidas no texto, é do maior interesse. 
Também a apresentação da recursão em dois exemplos, 
o primeiro — clássico — do cálculo do factorial de n, 
a o segundo da torre de Hanói, é de leitura simples e 
agradável. 


Smillie, Keith W, (1974). APL/360 With Statistical 
Examples. Addison-Wesley Publishing Company, Rea- 
ding, Massachusetts. 


Outra excelente publicação do mesmo autor, mais 
avançada em complexidade matemática e exigindo tam- 
bém mais noções de estatística, não muito profundas ou 
pelo menos superficiais, pois os conceitos são desenvol- 
vidos no texto, muitas vezes através de programas, 
como veremos. 

Uma vez que na crítica do outro manual do Profes- 
sor Smillia acabámos por falar do estudo que aí apre- 
senta da recursão, parece-nos interessante começar por 
referir aqui também o mesmo assunto. Este é agora 
abordado de maneira muito clara, com um capítulo com- 
pleto dedicado a esse tema. Se, por um lado, repete os 


135 


mesmos exemplos, outros surgem: é o caso dos números 
do Fibonacci, para os quais são apresentadas nada 
menos do que cinco funções geradoras: a primeira re- 
cursiva dupla, uma segunda chama outra recursiva sim- 
ples, a quarta interactiva, e a quinta, finalmente, gera 
todos os termos até uma ordem dada. Este capítulo de- 
dica-se ainda à função de Ackermann, ao cálculo dos 
dígitos de e numa base qualquer e a ordenações. 


No concernente ao estudo da estatística, as defini- 
ções surgem muitas vezes por via programática; de acor- 
do com a nossa própria experiência, isso é ideal para a 
experimentação com estes conceitos. Recordamos aqui, 
de passagem, a introdução das medidas estatísticas de 
localização (a média, a moda e, em especial neste con- 
texto, a mediana) e de dispersão. A moda só surge muito 
mais tarde do que é habitual, não obstante ser dos con- 
ceitos estatísticos de apreensão mais simples e natural, 
já depois do estudo das tabelas de contingência. Outros 
temas de estatística incluem histogramas, coeficientes 
de correlação simples e parciais, regressões lineares 
simples e múltiplas, o teste de Wilcoxon, a análise da 
variância e cálculo de probabilidades discretas. 


Estranhámos ver no decurso da exposição o autor 
usar uma função de inversão de matrizes escrita em APL, 
em vez da função nativa da linguagem na versão APL/ 
[360 usada, mas pensamos que isso terá sido motivado 
pela necessidade de apresentar um bom exemplo de uso 
das outras funções aplicadas a quadros. Também pode- 
riamos esperar num livro destes, datado do ano passado, 
o uso de certas extensões da linguagem (o APLSV foi já 
usado pelo autor noutras publicações); talvez isso tenha 
sido unicamente por uma limitação de espaço ou de 
tempo, mas sem dúvida melhoraria ainda mais a imagem 
do trabalho. 

Alguns capítulos são notáveis. Falámos já do da 
recursão. Outro é dedicado à simulação, de importância 
enorme no âmbito da estatística, e aí lemos com imenso 
agrado o tratamento do teorema do limite central e a sua 
comprovação empírica no caso da convergência da lei 
binomial para a normal, através de tiragens aleatórias 
simuladas pelo próprio APL, via função geradora de nú- 
meros aleatórios, Tal exposição é do maior valor peda- 
gógico e permite ver com clareza um resultado funda- 
mental em toda a teoria das probabilidades e da esta- 
tística. 


Walker, Terry (1972). An Introduction to Computer 
Science: APL Language Programming. Allyn and Bacon, 
Inc., Boston, Massachusetts. 


Um livro introdutório, talvez comparável com o de 
Prager (1971), incluido numa série, todos do mesmo 
autor, onde os outros livros tratam ou de outras lingua- 
gens ou de conceitos vários de computadores. 


Da facto pensamos dever esta obra ser lida em 
conjunção com os restantes livros de Walker, visto tratar- 
-se de uma colecção com carácter de referência, embora 
podendo ser também seguida para fins de ensino. No que 
respeita ao APL muitos tópicos poderiam ter sido apro- 
fundados, o que a verificar-se do mesmo modo nos res- 
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tantes livros tornaria, porém, esta colecção de certo 
modo enciclopédica. 


Weidmann, Clark (1975). Handbook of APL Progra- 
mming. Petrocelli/Charter, New York. 


Devido ao facto de este livro ser de publicação 
muito recente não tivémos até esta data possibilidade de 
obter uma cópia. Incluímo-lo aqui para tornar esta lista 
tão completa quanto possível (repetindo, o que dissémos 
no início deste artigo, não ser a mesma exaustiva), pois 
as críticas preliminares publicadas em revistas de infor- 
mática garantem-nos estar este livro de Weidmann pre- 
destinado a grande projecção no domínio do APL. 


Colloque APL. Institut de Recherche d'Informati- 
que et d'Automatique (IRIA), Domaine de Voluceau, Roc- 
quencourt, França (1971). 


Estas são as actas do primeiro congresso do APL 
realizado na Europa, em Setembro de 1971. Abrangem 
excelentes artigos, a par de outros menos bons, agrupa- 
dos em quatro capítulos: ensino e aplicação da lingua- 
gem APL (dois capítulos), outro sobre aspectos teóricos 
a o quarto acerca da implementação da linguagem, 

No primeiro capítulo podemos citar um artigo de 
Iverson, reproduzido em apêndice em lverson (1972). 
sobre o uso do APL como notação algébrica. Outros ar- 
tigos referem-se à utilização da linguagem no ensino da 
engenharia, da física, da estatística, da investigação ope- 
racional, etc. 

O capítulo sobre os aspectos teóricos contém três 
artigos, de grande profundidade qualquer deles. Um, de 
Abrams, tratando de aspectos formais da semântica onde, 
entre outras teses é abordada uma nova maneira de in- 
dexar quadros, baseando-se na redução por um operador 
novo, a introduzir, dos elementos a indexar; tal teoria, 
notável há quatro anos e, por exemplo, implementada no 
APL/1100, parece-nos vir a ser ultrapassada em estru- 
turas de quadros generalizados, O segundo artigo, sobre 
o mesmo tópico, é de Robinet (compare-se com o que 
dissemos acerca de Abrams et Lacourly (1972), Robinet 
(1971) e Perlis (1975), acima). Finalmente o terceiro 
artigo deste capitulo é de Lathwell e Mezei (outros dois 
implementadores), a descrição formal da linguagem em 
termos da própria linguagem, fundamental nara todos os 
estudos teóricos sobre os fundamentos do APL, e hoje 
largamente usada pelos implementadores de interpreta- 
dores e compiladores e construtores de máquinas APL. 

segundo sabemos este trabalho foi o primeiro onde 
foram usados sistematicamente os overadores de eixos 
(catenação — laminação, expansão, compressão, reflexão 
rotação, subida, descida, redução e varrimento), para os 
quais se prevê a definição de uma sintaxe uniforme a 
todos eles. É curioso que aqui o operador de varrimento 
(introduzido apenas em 1973 com o APLSV) só foi de- 
finido para a soma aritmética e igualdade e desigualdade 
lógicas; parecia óbvio tê-lo estendido a outras funções. 
No respeitante aos operadores de subida e descida ao 
longo dos eixos, para quadros, ainda não existem hoje 
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implementações práticas, segundo julgamos saber, em- 
bora haja rumores de eles terem feito parte do reportó- 
rio de instruções do STAR (compare-se com o que foi 
dito em Demars, Rault et Ruggiu (1974), acima). Este 
artigo usa extensivamente uma função de IF lógico, 
como aliás fazem outros autores (e. g. Smillie), imple- 
mentada com base na compressão pela expressão lógica 
e, por isso, de incorporação pouco provável na linguagem. 

Outro prato forte deste congresso foram, já então, 
as máquinas APL, a que são dedicados três artigos no 
último capítulo. 


APL Congress 1973. Gjerlov, P., Helms, H. J. and 
Nielsen, J., editors. North-Holland/American Elsevier, 
Amsterdam-London-New York (1973). 


Estas actas do segundo congresso europeu do APL 
contêm 65 artigos sobre variadíssimos tópicos, estudos 
teóricos a aplicações, desde a codificação em APL de 
trechos musicais à programação gráfica. Todos os artigos 
estão em Inglês. 

O congresso foi realizado sob o signo de algumas 
críticas violentas dirigidas à linguagem. Entre elas são 
notáveis as formuladas pelo Professor Dijksta na recep- 
ção do prémio Turing de 1972 (Dijkstra, E. (1972). «The 
humbls programmer». Comm. of the ACM, 15, Vol 10, 
pg. 859 et ff.). Estas críticas relacionam-se de perto com 
dois aspectos: o primeiro de repúdio à tendência de cer- 
tos autores para as linhas de código excessivamente 
compactas (as «one liners»), que por vezes abarcam 
todo um programa; compare-se com o que diz Ársac no 
prefácio a Robinet (1971). O segundo aspecto das cri- 
ticas de Dijkstra é sobre a falta de estruturas lógicas 
de controlo da execução das funções, dos tipos IF... 
THEN... ELSE, DO... END, e DO WHILE... END, e al- 
guns dos trabalhos do congresso de Copenhaga são de- 
dicados a possíveis saneamentos desta situação. 

Algumas semanas antes do congresso havia apare- 
cido porém um artigo de Falkoff e Iverson (Falkoff, A. 
and Iverson, K. E. (1973). «The design of APL». IBM J. 
of Research and Development, 17, pp 324-334), onde 
são precisamente criticadas essas estruturas de controlo 
por serem dependentes do contexto e estarem limitadas 
a escolhas binárias em cada passo. Entretanto foi anun- 
ciado também o APLSV e esta linha de investigação pa- 
rece ter sido, por agora pelo menos, relegada para plano 
secundário. Notamos aqui, de passagem, ser possível si- 
mular a instrução composta IF... THEN... ELSE muito 
facilmente por meio da função de execução do APLSV, 
mas essa utilização, no sentido da programação estrutu- 
rada, transparece como estranha à estrutura do APL. 

Outra linha de desenvolvimento futuro hoje muito 
em destaque (veja-se acima a crítica a Demars, Rault et 
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Ruggiu (1974)) diz respeito aos quadros generalizados, 
funções e operadores, a que nos referimos múltiplas 
vezes; também nestas actas foram dedicadas várias co- 
municações a estes assuntos. Pensamos que esta orien- 
tação será das mais importantes e trará, de futuro, as 
modificações mais drásticas a introduzir no APL. Para 
já permitirá a existência de funções definidas multiádi- 
cas, sendo hoje uma das limitações da linguagem a va- 
lência das funções só poder ser O (para as definidas) 
o 1 ou 2 (para as nativas e para as definidas), outra 
séria limitação, aliás notada por Robinet (1971) e De- 
mars et al, (1974). 


Ainda uma terceira linha de investigação orienta-se 
no sentido de processamentos com pares de valores, su- 
jeitos no cálculo às leais das operações com números 
complexos, e a ela é dedicado um artigo; é curioso ter 
este aspecto sido abordado também nesse artigo de Fal- 
kofí e Iverson, mas sem ser sugerida saída alguma. So- 
mos dos que pensam só vir este tópico a ser solucionado 
num âmbito mais geral, eventualmente após terem sido 
implementados os quadros de quadros; outras lingua- 
gens, como o PL/I|, o ALGOL 68 e o próprio FORTRAN 


tratam-no sem grandes complexidades, é certo, mas sem 
qualquer generalidade de interessa. 


APL 75: Congress Proceedings. Association for Compu- 
ting Machinery, New York (1975). 


Meia centena de artigos é o activo do terceiro con- 
gresso APL na Europa, subordinado aos temas: desen- 
volvimento, implementação e aplicações da linguagem. 

Uma boa parte deles foi novamente dedicada à in- 
terpretação directa pelo «hardware». Alguns outros tópi- 
cos de interesse e sem dúvida de grande relevo futuro 
dizem respeito à verificação dos programas (tema já 
abordado também em 1973 e noutros artigos dispersos 
pela literatura e da importância quanto à própria fiabili- 
dade das máquinas APL), redes de computadores e mon- 
tagem 6a acesso a bases de dados em APL. Ainda outra 
tendência, já verificada, no sentido de criticas acerbas 
a certas estruturas e implementações da linguagem vol- 
taram à ordem do dia. 


Parece, na realidade estar o APL a despertar um 
interesse de uma parte da comunidade informática mun- 
dial que ainda há escassos anos seria um sonho imagi- 
nar. Em particular é agora urgente uma normalização da 
linguagem; se a sua falta permitiu uma grande flexibili- 
dade aos implementadores durante vários anos, possibi- 
litando-lhes delinear as melhores opções antes de as tor- 
nar disponíveis, a proliferação de sistemas e, daqui em 
diante também de máquinas, torna-a indispensável, como 
condição para a expansão continuada e harmónica do 
APL. 
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Resumo dos Artigos publicados na «Técnica» n.º 43] 
Ano L — Dezembro 1975 


C. D. U. 621.373:621.316.825.4 


J. F. BORGES DA SILVA 
J. M. VIEGAS GONÇALVES 


OSCILAÇÕES AUTO-EXCITADAS NUM 
CONDENSADOR-TERMISTOR. 


Técnica N.º 431 — L— 12-1975, p. 89 a 102 


CIRCUITO 


Faz-so a análise do comportamento dinâmico da associação 
em paralelo dum condensador com um termistor alimentada 
por um gerador ideal de corrente, e em particular da possi- 
bilidade de ocorrência de oscilações auto-excitadas. Compa- 
ram-se resultados experimentais com os resultados da 
simulação numérica do modelo teórico proposto. São des- 
critas os experiências que permitem determinar os perâme- 
tros do termistor a uma montagem destinada a simular os 
elevados valores de capacidade que foram necessários nas 
experiências. Utilizando a simulação numérica, faz-se o 
estudo o interpretação fisica do comportamento do sistema 
para diferentes valores da capacidado do condensador. 
Apresenta-se, para caracterização do termistor no dominio 
incremental, um esquema eléctrico equivalente deduzido das 
equações diferenciais não-linearas que regem O seu com- 
portamento, 


C. D. U. 624.07:518.3 
E. R. DE ARANTES E OLIVEIRA 


REFLEXÕES SOBRE A EVOLUÇÃO DA ENGENHARIA 
DE ESTRUTURAS. 


Técnica No 431 — L— 12-1975, p. 103 a 116 


Recordam-sa as duas revoluções sofridas pela Engenharia 
de l|Estruturas, provocadas respectivamente nela utilização 
dos princípios da Mecânica, no século XVIII, e dos compu- 
tadores digitais em meados do século XA. O desenvolvi- 
mento da técnicas de discretização, como a dos elementos 
finitos, velo permitir a resolução de numerosos problemas. 
Sente-se no entanto a necessidade de uma profunda modifi- 
cação dos métodos tradicionais que torne possivel uma mais 
eficaz exploração das potencialidades das máquinas. 
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€. D. U. 621.374 


L. T. MAGALHÃES 
M. P. O. RICOU 


MAQUINA DE TURING PARA CONVERTER EXPRESSÕES 
EM NOTAÇÃO “INFIX” EM EXPRESSÕES EM NOTAÇÃO 
"POSTFIX" — |. DEFINIÇÃO DA MAQUINA E REDUÇÃO 
DE ESTADOS. 


Esto artigo é o primeiro de dois dedicados ao proiecto de 
uma máquina de Turing para converter notação intix em 
postfix. São amplamente utilizados, como auxiliares no 
prójecto, métodos de análise de estrutura algébrica de maá- 
quinas sequenciais, que têm vindo a ser estabelecidos no 
dominio da Teoria dos Automata. 

Nesto primeiro artigo é apresentado um algoritmo, 
definida uma máquina de Turing nele baseada, e reduzidos 
os estados da máquina de modo a obter uma máquina minima. 


J. MARQUES HENRIQUES 
UMA BIBLIOGRAFIA DE APL, 
Técnica N.º 431 — L— 12-1975, p. 127 a 138 


Esto artigo é uma colectânea de breves críticas bibliográfi- 
cas de livros sobre a linguagem APL, À variedade de tópicos 
tratados nestas obras reflecte a profunda importância assu- 
mida por esta linguagem em todos os domínios da informá- 
tica. 
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Synopsis of articles published in «Técnica» n.º 431 
L —Dezember 1975 


Cc. D. U. 621.374 


L, T. MAGALHÃES 
M. P. O. RICOU 


TURING MACHINE FOR THE TRANSLATION FROM INFIX 
TO POSTFIX NOTATION — |. DEFENITION AND STATE 
REDUCTION. 


Técnica N.º 431 — L-— 12.975, p. 117 a 126 


This is the first of two papers related to the design of a 
Turing machine for the translation of infix notation to postfix 
notation. Methods for the analysis of algebraic structure of 
sequential machines, recently developed as part of Automata 
Theory, are largely used. 

im this first paper we present an algorithm, define a 
Turing machine realizing it, and reduce machine states in 
order to get a minimal machine, 


C. D. U. 681.14.016 
J. MARQUES HENRIQUES 
A BIBLIOGRAPHY OF APL. 
Técnica N.º 431 — L— 12-1975, p. 127 a 138 
This survey paper brings togather a collection of brief ravi- 
ews of books on APL. The manifold of topics dealt with 


im theso works shows tha primary importance achieved by 
the language in all fialds of information processing. 
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C. D. U. 621.373:621.316.825.4 


J. F. BORGES DA SILVA 
J. M. VIEGAS GONÇALVES 


SELF-EXCITED OSCILLATIONS IN A CONDENSER-THER 
MISTOR CIRCUIT. 


Técnica N.º 431 — L— 12-1975, p. 89 a 102 


The dynamic behaviour of a condenser and a thermistor 
connected in parallel with an ideal current souce is analysed, 
particular attention being given to the possibility of self- 
-excited oscillarions. Results from experiments are compared 
with results obtained by numerical simulation of the theo- 
retical model proposed. Experiments required for the deter- 
mination of thermisltor parameters and a method for the 
simulation of large capacitance values are described. System 
behaviour was further investigated for various capacitance 
values using numerical simulation, tha results being interpre- 
ted physically. à thermistor small signal equivalent circuit, 
derived from the non-linear differential equations describing 
its behaviour is presented. 


C. D. U. 624.07:518.3 


E. R. DE ARANTES E OLIVEIRA 


CONSIDERATIONS ON THE EVOLUTION OF STRUCTURAL 
ENGINEERING. 


Técnica N.º 431 — L-— 12-1975, p. 103 a 116 


Tha two revolutions suffered by Structural Engineering, res- 
pectively in tha XVIIL th century, due to the introduction 
ot tha principles of Mechanics, and in the XA th century, 
due to tha advent of computers, are first remembered. 
Attention is then called to tha fact that a fully efficiant 
exploration of the capacity of computers can be possible 
only if à deeper change of traditional methods is achieved. 


Artigo N.º Data Pág. 
6 —. CIÊNCIAS APLICADAS 
62 — ENGENHARIA 


620 —. ENSAIOS E ACÇÕES FÍSICAS SOBRE OS MA- 
TERIAIS 


C. D. U. 620.193:66.046.37:621:438:546.31.9 


«Heigh — Temperature corrosion 
of nickel-base alloys in the 
gas turbine» ........ psciicaaa d 427 6-975 317 


César Augusto Correia de Sequeira 


621 — ENGENHARIA MECÂNICA 


C. D. U. 621:438:546.31.9:620.193:66.046.37 


«Heigh — Temperature corrosion 
of nickel-base alloys in the 


gas turbine» ..... a ai a 427 6-975 317 


César Augusto Correia de Sequeira 


621.35 — ENGENHARIA ELECTROTÉCNICA 


C. D. U. 621.3:517.512.2 


Estudo experimental da aplicação 
do cepst.o (1) à detecção de 
ecos e periodicidades em si- 


Gi RR IRES 424 3-975 167 


Luis Henrique Martins Borges de Almeida 
Guilherme Diniz Moreno da Silva Arroz 


C. D. U. 621.3:517.512.2 
«Cepstrum» 
Estudo experimental da sua apli- 
cação ao sinal falado «Voco- 


DOME. sucursais id 428 7-975 389 


lriarte José Araújo Esteves 
José Alberto Teixeira 


C. D. U. 621.317.34 
Identificação de sistemas lineares 


po: amostragem finita no do- 
mínio do tempo ............... 423 2-975 133 
J. M. Viegas Gonçalves 

C. D. U, 621.374 


A utilização de variáveis pseudo- 
-aleatórias na identificação de 
ERSRUCREARSA e qicicie a 427 6-975 343 

S. D, Antunes 

C. D. U. 621.392.1:621.396.029.1 


Analysis of fade durations on dif- 
ferent microwave paths ...... 421 


José A. C. Saraiva Mendes 


10-974 23 
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C. D. U. 621.396.65 


Estudo da propagação da banda 
dos 11 GHz para ligações por 
Doo SN SOR PR 426 


José A. C. Saraiva Mendes 


5-975 289 


621.9 — FERRAMENTAS. MÁQUINAS-FERRAMENTAS. 
OPERAÇÕES MECÂNICAS 


C. D. U. 621.9:338.45 


Produção de máquinas-Ferramen- 


tas em pequena série ...... 422 1Y1-974 89 


Luciano de Oliveira Faria 


622 —. EXPLORAÇÃO MINAS E TRATAMENTOS Ml- 
NEIROS 


C. D. U. 622.003 
Bases do estudo económico de um 


empreendimento mineiro ... 423  2-975 111 


C. Dinis da Gama 


624 — ENGENHARIA CIVIL 


C. D. U. 624.022.327.002.22:699.841 


Metodologia para o estudo sísmico 
de edifícios prefabricados por 
EOMRRCNDES serias DOS ira 423 
S. Ferry Borges 
A. Saraiva 
J. Pereira 
V. Monteiro 


2-975 121 


C. D. U. 624.025.3.699.844 


Pavimentos flutuantes. Ensaio em 
modelo para qualificação do 
isolamento à transmissão de 


ruídos de percussão ......... 427 6-975 349 


P. Martins da Silva 
C. D, U. 624.04:518.6 
O método dos gradiantes conjuga- 


dos e a sua aplicação na aná- 


lise de estruturas ........c.. 421 10-974 1 


E. R. de Arantes e Oliveira 
L. F. Gonçalves de Azevedo 
C. D. U, 624.072.23 


Vigas continuas parcialmente en- 
castradas ........ is VA 


Jd. S. Brazão Farinha 


3-975 185 


C. D. U. 624.074.6 


Sobre torres espiadas —- Método 
de cálculo e um exemplo de 
Apiicação: icssssssescicsaiacos 422 11-974 49 


Luis Fernando Gonçalves Morgado de Azevedo 


Artigo No Data Pãg. 


C. D. U. 624.131.388:627.82 


Ensaios de fundações rochosas de 
barragens e seu comporta- 


PORTÃO: mussi aaa ie da dc Ada 428 7-975 369 


J Laginha Serafim 
Marciano L. Guerreiro 


627.8 — BARRAGENS 


C. D. U, 627.82:624.131.388 


Ensaios de fundações rochosas de 
barragens e seu comporta- 


ici gds (QI 428 1-975 369 


J, Laginha Serafim 
Marciano L. Guerreiro 


628 — ENGENHARIA E TECNOLOGIA SANITÁRIAS 


C. D. U. 628.5:551.510.4 


Um modelo gaussiano para a dis- 
persão de poluentes na at- 


Seis ici qo PRDC q PDR 428 7-975 383 


Jd. J. Delgado Domingos 
M. Teresa Pontes 


C. D. U. 628.52:533.15:551.,554 


Diffusion dans les couches limites 


atratifiãos .:.iiiccisisasecsessei 424 3-975 221 


8 
425  4-975 237 


Jd. Mathieu 


629 —. CONSTRUÇÃO DE VEÍCULOS — ASTRONAU- 
TICA 
C. D. U. 629.19:001+62 


Spaco and scientific and techno- 


logical progress ............... 426 5.975 305 


O. M. Belotserkovskii 


C. D. U. 629.19:681.14 


Tratamento, num computado: IBM 
1130, da informação acerca 
da trajectória do um satélite 
artificial ....... Par SERRA TIS 422 11-974 97 


Teresa Levy e Almeida Loureiro 
C. D. U. 629.19:551.501 


Obtenção dos dados meteorológi- 
cos através de um satélite ar- 
REG quieta rondis ads 421 10-974 33 


Bettencourt Faria 


Artigo N. Data Pág. 


C. D. U. 629.19:681.14:551.501 
Conversor digital análogo temporal 423 2-975 141 
Rodrigo Martins 
Isaura Martins 
66 — TECNOLOGIA QUÍMICA 
C. D. U. 66.046.37:621:438:546.31.9:620.193 


«Heigh — Temperatura corrosion 
of nickel-base alloys in the 
gas turbinasS cseesssanemees 427 6-975 317 


César Augusto Correia de Sequeira 
C. D. U. 663.63 


Reutilização de água .........cc. 427 6-975 355 


J. Sarafana 


681.14 — MÁQUINAS DE CALCULAR 
C. D. U. 681.14:517.5(084) 


Um programa para o traçado de 
curvas de nivel em pequenos 


calculadoros saci. 428 7-9756 409 


J. !, Amarante dos Santos 
C. D. U. 681.14:551.501:629.19 
Conversor digital análogo temporal 423 2-975 1419 
Rodrigo Martins 
Isaura Martins 
C. D. U. 681.14:629.19 
T.atamento, num computador IBM 
113€, da informação acerca 


da trajectória de um satélite 
MENERA ea ap 422 11-974 97 


Teresa Levy 
Almeida Loureiro 


69 —. CONSTRUÇÃO E EDIFÍCIOS, MATERIAIS E PRO- 
CESSOS DE CONSTRUÇÃO 


699.8 —. PROTECÇÃO DAS CONSTRUÇÕES 
C. D. U. 699.841:624.022.327.002.22 


Metodologia para o estudo sísmico 
de edifícios prefabricados por 
Cio fa E PRP VA PESE 423 2-975 121 
- Ferry Borges 
- Saraiva 
Pereira 
- Monteiro 


xt 


C. D. U. 699.844:624.025.3 


Pavimentos flutuantes. Ensaio em 
modelo para qualificação do 
isolamento à transmissão de 
ruídos de percussão ......... 427 


P. Martins da Silva 


6-975 349 
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INSTITUTO SUPERIOR TÉCNICO 


ÍNDICE POR ASSUNTOS DOS ARTIGOS PUBLICADOS 


NUMEROS 421 A 428 VOLUME XXXVI 


Artigo N.º Data Pãg. 


OC — FUNDAMENTOS GERAIS DO CONHECIMENTO 
E CULTURA 


001 — CIÊNCIA E CONHECIMENTO EM GERAL. SA- 
BER. ERUDIÇÃO 


Cc. D. U. 001+62:629.19 


Space and scientific and technolo- 


gical Prograss ..amccccmssess 426 5-975 305 


O. M. Belotserkovskii 


007 — ACTIVIDADE E ORGANIZAÇÃO EM GERAL: 
INFORMAÇÃO, COMUNICAÇÃO E TEORIA DO 
CONTROLE EM GERAL, CIBERNÉTICA, INCLUIN- 
DO «ENGENHARIA HUMANA» 


C. D. U. 007 


Memo-ando sobre a inteligência 


ErTIGIO] cocarecsecarinesseaníria MT 6-975 329 


Helder Coelho 


3 —. CIÊNCIAS SOCIAIS 


338 —. PRODUÇÃO. CRIAÇÃO DA RIQUEZA. SITUA- 
ÇÃO ECONÔMICA 


C. D. U. 338.45:621.9 


Produção de máquinas-ferramentas 


em pequena série ........ 422 11-974  B9 


Luciano de Oliveira Faria 


43 — LÍNGUAS GERMÁÂNICAS 


C. D. U. 438:546.31.9:620.193:66.046.37:621 


«Heigh — Temperature corrosion 

of nickel-base alloys in the 

gas tuiDInAM susssnessssreãs 427 6-975 317 
César Augusto Correia de Sequeira 
5 —. CIÊNCIAS PURAS 


5.1 — MATEMÁTICA 
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OUTUBRO DE 1974 A JULHO DE 1975 


Ártigo N.o Data Pão. 


C. D. U. 517:537.313 


Sistemas de segunda ordem — 
Resposta no tempo e na fre- 
quência ............ curar MB 4-975 277 


Ricardo Charters d'Azevedo 


C. D. U. 517.5(084):681,14 


Um programa para o traçado de 
curvas de nível em pequenos 


calculadores ......cciiicesessas 428 7-975 409 


J !. Amarante dos Santos 


C. D. U. 518.6:624.04 


O método dos gradiantes conju- 
gados e a sua aplicação na 
análise de estruturas ......... 421 10-974 1 

E. R. de Arantes e Oliveira 


L. F. Gonçalves de Azevedo 


C. D. U. 518.6:539.3) 


Um método de forças ficticias 
para a resolução de proble- 
mas geometricamente não li- 
DONE area a 421 
E R. Arantes e Oliveira 
R. Teixeira Duarte 


19-974 T 


C. D. U. 517.512.2:621.3 


Estudo experimental da sua apli- 
cação ao sinal falado. «Voco- 
DOE assis aus an TA 428 
lriarte José Araújo Esteves 
José Alberto Teixeira 


7-975 389 


C. D. U. 517.512.2:621.3 
Estudo experimental da aplicação 
do Cepstro (1) à detecção 
de ecos e periodicidades em 
sinais ......... araiçédaa o. 4248 
Luis Henrique Martins Borges de Almeida 
Guilherme Diniz Moreno da Silva Arroz 


C. D. U. 518.6 


3-975 167 


Solução da equação de Poisson 
pelo método de ritz-processo 


geral para qualque: contorno 423 2-975 127 


Álvaro Bento Leal 


Artigo N.º Data Pãg. 
53 — FISICA 


C. D. U. 531.7.083:533.6.015 


Measurements in the region of 


recirculation behind a disc... 426 5-975 297 


D. F. G. Durão 
Jd. H. Whitelaw 


C. D. U. 532.5:599.3 


A unified theory of small distur- 
bances within compressibla 
lamina: and tu-bulent boun- 
dary layers .ecciiemeseeeciteso 422 11-974 19 


G. R. Inger 


C. D. U, 533.15:551.554:628.52 


Diffusion dans les couches limites 
SUStHTIÕOS icacissiscmiiaa casi: 424 3-975 221 
e 


425 4-975 237 
Jd. Mathieu 


C. D. U. 533.6.015:531.7.083 
Measurements in the region of re- 
circulation behind a disc ... 426 
F. G. Durão 
H. Whitelaw 


5-975 297 
D. 
Jd. 


C. D. U. 537.313:517 


Sistemas de segunda ordem — 
Respostas no tempo e na fre- 
QUERIA Capssias se isa ra inca Mni rod dai 425  4-975 277 


Ricardo Charters d'Azevedo 


C. D. U. 539.232.543.42 


Espectros de massa e energia na 
formação de películas finas 
usando o processo de «eva- 
poração por impacto» ....... 424 3975 213 


Leopoldo J. M. Guimarães 


C. D. U. 539.3:532.5 


A unified theory of small distur- 
bances within compressible 
laminar and turbulent boun- 
dáty EIGIE asassamesmecenars 422 11-974 79 


G. R. Inger 
C. D. U. 539.31:518.6 


Um método de forças fictícias 
para a resolução de proble- 
mas geometricamente não li- 
TIOATOS  usunssenisiniiscascssõo 421 10-974 7 


E. R. Arantes e Oliveira 
R. Teixeira Duarte 


àrtigo N. Data Pão. 


C. D. U. 539.56 
Introdução à mecânica da fractura 421 10-974 13 


Salgado Prata 


54 — QUÍMICA-MINERALOGIA 
546 — QUÍMICA INORGÂNICA 


C. D. U. 546.31.9:620.193:66.046.37:621:438 


«Heigh — Temperature corrosion 
of nickel-base alloys in the 
gas turbinas .....cceseasesesess 427 6-975 317 


César Augusto Correia de Sequeira 


55 — GEOLOGIA E CIÊNCIAS AFINS — METEOROLO- 
GIA 


C. D, U. 551.601:629.19 


Obtenção dos dados meteorológi- 
cos através de um satélite a:- 
CERDEIE pESanAsa SE EMA RETA A AT ÇÃ 421 10-974 33 


Bettencourt Faria 


C. D. U. 551.501:629.19:681.14 
Conversor digital análogo temporal 423 2-975 141 


Rodrigo Martins 
Isaura Martins 


C. D. U. 551.510.4:628.5 


Um modelo gaussiano para a dis- 
persão de poluentes na at- 


A ca a A TE Sad 428 7-975 383 


J. J. Delgado Domingos 
M. Teresa Pontes 


C. D. U. 551.554:628.52:593.15 


Diffusion dans les couches limites 
GtratHIGOS usessiassarsa rasas 424 3-975 221 
8 
425 4-975 237 
J. Mathieu 


C. D. U. 651.577:551,094.7 


Analysis of rain rate distributions 
in Portugal and its applica- 
tion to the microwave propa- 
gation above 10 GHZ ...... 425 4-975 263 


José A. Saraiva Mendes 
Carlos Monteiro 
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— BOLETIM — COLLECTION 
Sociedade Geológica de Portugal Czechoslovak Academy of Sciences 
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— DYNA 


Instituto Nacional de Racionalizacion y Normalizacion 


— BOLETIM DA ELECTROQUIMICA E CORROSÃO | 
Instituto Nacional de Investigação Industrial — EL ESTANO Y SUS APLICACIONES 


Tin Research Institute 


Asociacion de Engenieros Industriales de Espafia 


— BOLETIM 
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— ECONOMIC REVIEW 
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— FUNDEXPORT 
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Revista da Sociedade Portuguesa de Geotécnica 
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— INDÚSTRIA DE MOÇAMBIQUE 
Associação Industrial de Moçambique 
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— METALLURGIA AND METAL FORMING 
National Association of drop Forgers and Stampers 


— MACHINE TOOL REVIEW 


— NOTICIÁRIO DOS ELASTÔMEROS 
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— ONU (CHRONIQUE MENSUELLE) 
Servico de L'information de L'organisation des Na- 
tions Unies 
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Câmara de Comércio e Indústria Luso-Alemã 
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Office National D'Études et de Recherchse Aéros- 
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— REVISTA DE METALURGIA 
Cenim - Centro Nacional de Investigaciones Meta- 
lurgicas 


— RAIL INTERNATIONAL 
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de fer et de L'Union Internationale des Chemins 


do fer 


— REVISTA TÉCNICA SULZER 


— REVISTA DANFOSS 


— REVISTA AGRONÔMICA 
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— TECHNIQUES FRANÇAISES (BATIMENT TRAVAUX 
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— TECHNIQUES FRANÇAISES (INDUSTRIE DE 
METAUX 
Services Economiques de L'Embassade de France 
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— TECHNIQUES FRANÇAISES (CHIMIE GAZ PETROLE) 
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ERRATA TÉCNICA 430 


Termodinâmica molecular do gás perfeito 


67 — Na penúltima linha do «Abstract», em vez 
de exprerimental deve ler-se experimental. 


— Na linha 11 da Introdução, onde está 
termodinêmicas deve ler-se termodiná- 
micas. 


ãg. 68 —- Na linha 4 depois da equação (5) deve 


ler-se e, em vez de E. 


— À equação (8) deve ler-se 
w,e — e;/kT 
Z 


— Na linha 5 após a figura 1 deve ler-se 
de, : 


70 —- Na linha 3 depois da expressão (27), 
onde está N, deve ler-se vi 


— O 2,º membro da equação (28) é um pro- 
duto iterado. 


— Na equação (30) falta um parêntesis 
recto imediatamente a seguir ao último 
parêntesis curso, 


71 — Na linha 4, onde está B, deve lar-se B, : 


—— Na linha 2 depois da equação (47), onde 
v deve ler-se py. 


72 — Na expressão (48) não deve constar o 
último parêntesis curvo. 


— (O numerador do 2.º membro da equação 
(53) deve ler-se 


VB6 h 


74 — Na linha 4 da expressão (82), onde está 
'p deve ler-se “pj : 


— Na linha 3 da expressão (83), logo após 
o parêntesis recto inicial, em vez de P; 
deve ler-se “q, . 


J. A. MARTINHO SIMÕES e 
C. A. NIETO DE CASTRO 
Laboratório de Quimica-Fisica 
Instituto Superior Técnico 


— Pág. 75—0 índice do 2.º somatório da equação (85) 
é i<j<k. 


— Na linha 8 da quação (85), antes do pa- 
rêntesis recto, deve ler-se 


m+]] 


— A expressão de “p das equações (88) 
deve ler-se 


'p=4 98 el... 
— Na linha 2 da nota (6) deve ler-se 
a 
Cj=aTico 
— Pág. 76 —- Na equação (96) deve ler-se 
b à a 
ra dr; 9; 


— Na equação (97) deve ler-se 
per a 
rd CD 


— Na penúltima linha do parágrafo 3 (Cál- 
culos) deve ler-se C,em vez de Co- 


— Pág. 78 — Na tabela Il, linha 2, coluna 5, deve ler-se 
13.5 x 1074 


—— Na nota (c), onde está [5] deve ler-se 
[15]. 


— Na tabela |, linha 1, coluna 7, deve ler-se 


B 
&; 


— Pág. 79 — Na tabela III, linha 1, coluna 3, deve ler- 


se |. 
E + (C, | INk 


LOCO nte aa 0 
0000000000 0000000000: | 
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